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In primary school, we were told that there are four phases of matter: solid, liquid, gas, and 
plasma. In college, we learned that there are much more than four phases of matter, such as 
hundreds of crystal phases, liquid crystal phases, ferromagnet, anti-ferromagnet, superfluid, etc . 
Those phases of matter are so rich, it is amazing that they can be understood systematically by the 
symmetry breaking theory of Landau. However, there are even more interesting phases of matter 
that are beyond Landau symmetry breaking theory. In this paper, we review new “topological” 
phenomena, such as topological degeneracy, that reveal the existence of those new zero-temperature 
phases - topologically ordered phases. Microscopically, topologically orders are originated from the 
patterns of long-range entanglement in the ground states. As a truly new type of order and a truly 
new kind of phenomena, topological order and long-range entanglement require a new language 
and a new mathematical framework, such as unitary fusion category and modular tensor category 
to describe them. In this paper, we will describe a simple mathematical framework based on 
measurable quantities of topological orders {S, T, c) proposed around 1989. The framework allows 
us to systematically describe all 2-l-lD bosonic topological orders {i.e. topological orders in local 
bosonic/spin/qubit systems). 
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I. INTRODUCTION 

A. Orders, phase transitions, and symmetries 

Condensed matter physics is a branch of science that 
study various properties of all kinds of materials, such 
as mechanical properties, hydrodynamic properties, elec¬ 
tric properties, magnetic properties, optical properties, 
thermal properties, etc . Since there are so many dif¬ 
ferent kinds of materials with vastly varying properties, 
not surprisingly, condensed matter physics is a very rich 
field. Usually for each kind of material, we need a differ¬ 
ent theory (or model) to explain its properties. So there 
are many different theories and models to explain various 
properties of different materials. 

However, after seeing many different type of theo¬ 
ries/models for condensed matter systems, a common 
theme among those theories start to emerge. The com¬ 
mon theme is the principle of emergence, which states 
that the properties of a material are mainly determined 
by how particles are organized in the material. Differ¬ 
ent organizations of particles lead to different materials 
and/or different phases of matter, which in turn leads to 
different properties of materials. 

Typically, one may think that the properties of a ma¬ 
terial should be determined by the components that form 
the material. However, this simple intuition is incorrect, 
since all the materials are made of same three compo¬ 
nents: electrons, protons and neutrons. So we cannot 
use the richness of the components to understand the 
richness of the materials. The various properties of dif¬ 
ferent materials originate from various ways in which the 
particles are organized. The organizations of the parti¬ 
cles are called orders. The orders (the organizations of 
particles) determine the physics properties of a material. 

Therefore, according to the principle of emergence, the 
key to understand a material is to understand how elec¬ 
trons, protons and neutrons are organized in the mate¬ 
rial. However, to develop a theory for all possible or¬ 


ganizations of particles, we need to have a more precise 
description/definition of organizations of particles. 

First, we need to find a way to determine if two orga¬ 
nizations of particles should be regard as the same (or 
more precisely, belong to the same class or belong to 
the same phase) or not. Here, we need to rely on the 
phenomena of phase transition. If we can deform the 
system (such as changing temperature, magnetic field, 
or other parameters of the system) in such a way that 
the state of the system before the deformation and the 
state of the system after the deformation are smoothly 
connected without any phase transition,then we say 
the two states before and after the deformation belong to 
the same phase and the particles in the two states are re¬ 
garded to have the same organization. If there is no way 
to deform the system to connect two states in a smooth 
way, then, the two states belong to two different phases 
and the particles in the two states are regarded to have 
to two different organizations. 

We note that our definition of organizations is a def¬ 
inition of an equivalent class. Two states that can be 
connected without a phase transition are defined to be 
equivalent. The equivalent class defined in this way is 
called the universality class. Two states with different 
organizations can also be said to belong to different uni¬ 
versality classes. We introduce a formal name “order” to 
refer to the “organization” defined above. 

Based on a deep insight into phase and phase transi¬ 
tion, Landau developed a general theory of orders as well 
as transitions between different phases of matter^^^. Lan¬ 
dau points out that the reason that different phases (or 
orders) are different is because they have different sym¬ 
metries. A phase transition is simply a transition that 
changes the symmetry. Introducing order parameters 
that transform non-trivially under the symmetry trans¬ 
formations, Ginzburg and Landau developed Ginzburg- 
Landau theory, which became the standard theory for 
phase and phase transition.^ For example, in Ginzburg- 
Landau theory, the order parameter can be used to char¬ 
acterize different symmetry breaking phase: if the order 
parameter is zero, then we are in a symmetric phase; if 
the order parameter is non-zero, then we are in a sym¬ 
metry break phase. The symmetry breaking phase tran¬ 
sition is the process in which the order parameter change 
from zero to non-zero. 

Landau’s theory is very successful. Using Landau’s 
theory and the related group theory for symmetries, we 
can classify all of the 230 different kinds of crystals that 
can exist in three dimensions. By determining how sym¬ 
metry changes across a continuous phase transition, we 
can obtain the critical properties of the phase transi¬ 
tion. The symmetry breaking also provides the origin of 
many gapless excitations, such as phonons, spin waves, 
etc., which determine the low-energy properties of many 
systems.^’® Many of the properties of those excitations, 
including their gaplessness, are directly determined by 
the symmetry. 

As Landau’s symmetry-breaking theory has such a 
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broad and fundamental impact on our understanding of 
matter, it became a corner-stone of condensed matter 
theory. The picture painted by Landau’s theory is so 
satisfactory that one starts to have a feeling that we un¬ 
derstand, at least in principle, all kinds of orders that 
matter can have. One starts to have a feeling of seeing 
the beginning of the end of the condensed matter theory. 

However, through the researches in last 25 years, a dif¬ 
ferent picture starts to emerge. It appears that what 
we have seen is just the end of beginning. There is a 
whole new world ahead of us waiting to be explored. 
A peek into the new world is offered by the discovery 
of fractional quantum Hall (FQH) effect.® Another peek 
is offered by the discovery of high superconductors.^ 
Both phenomena are completely beyond the paradigm of 
Landau’s symmetry breaking theory. Rapid and excit¬ 
ing developments in FQH effect and in high Tc super¬ 
conductivity resulted in many new ideas and new con¬ 
cepts. Looking back at those new developments, it be¬ 
comes more and more clear that, in last 25 years, we 
were actually witnessing an emergence of a new theme in 
condensed matter physics. The new theme is associated 
with new kinds of orders, new states of matter and new 
class of materials beyond Landau’s symmetry breaking 
theory. This is an exciting time for condensed matter 
physics. The new paradigm may even have an impact in 
our understanding of fundamental questions of nature - 
the emergence of elementary particles and the four fun¬ 
damental interactions.®^^® 


B. The discovery of topological order 

After the discovery of high T^. superconductors in 
1986,^ some theorists believed that quantum spin 
liquids play a key role in understanding high 
superconductors^^ and started to construct and study 
various spin liquids.Despite the success of Landau 
symmetry-breaking theory in describing all kind of states, 
the theory cannot explain and does not even allow the 
existence of spin liquids. This leads many theorists to 
doubt the very existence of spin liquids. In 1987, a spe¬ 
cial kind of spin liquids - chiral spin state^®’^® - was 
introduced in an attempt to explain high temperature 
superconductivity. In contrast to many other proposed 
spin liquids at that time, the chiral spin liquid was shown 
to correspond to a stable zero-temperature phase and is 
more likely to exist.®®® At first, not believing Landau 
symmetry-breaking theory fails to describe spin liquids, 
people still wanted to use the symmetry breaking theory 
to characterize the chiral spin state. They identified the 
chiral spin state as a state that breaks the time rever¬ 
sal and parity symmetries, but not the spin rotation and 
translation symmetries.^® However, it was quickly real¬ 
ized that there are many different chiral spin states (with 
different spinon statistics and spin Hall conductances) 
that have exactly the same symmetry, so symmetry alone 
is not enough to characterize different chiral spin states. 


This means that the chiral spin states contain a new kind 
of order that is beyond symmetry description.^^ This new 
kind of order was named^® topological order. ®®^ 

But experiments soon indicated that high-temperature 
superconductors do not break the time reversal and par¬ 
ity symmetries and chiral spin states do not describe 
high-temperature superconductors.^® Thus the concept 
of topological order became a concept with no experi¬ 
mental realization. 

Although the concept of topological order is introduced 
in a theoretical study, about a state that is not known to 
exist in nature, this does not prevent topological order 
to become a useful concept. As we will see later that the 
concept of topological order contains inherent self con¬ 
sistency and stability. If we believe in nature’s richness, 
all nice concepts should be realized one way or another. 
The concept of topological order is not an exception. 

Long before the discovery of high superconductors, 
Tsui, Stormer, and Gossard discovered FQH effect,® such 
as the filling fraction v = 1/m Laughlin state^® 

^ ( 1 ) 

where Zi = Xi + iyi. People realized that the FQH states 
are new states of matter. However, influenced by the 
previous success of Landau’s symmetry breaking theory, 
people still want to use order parameters and long range 
correlations to describe the FQH states. ^®^^® But, if we 
concentrate on physical measurable quantities, we will 
see that all those different FQH states have exactly the 
same symmetry and conclude that we cannot use Landau 
symmetry-breaking theory and symmetry breaking order 
parameters to describe different orders in FQH states. 
So the order parameters and long range correlations of 
local operators are not the correct way to describe the 
internal structures of FQH states. In fact, just like chiral 
spin states, FQH states also contain new kind of orders 
beyond Landau’s symmetry breaking theory. Different 
FQH states are also described by different topological 
orders.^® Thus the concept of topological order does have 
experimental realizations in FQH systems. 

In addition to the Laughlin states, more exotic non- 
abelian FQH states were proposed in 1991 by two inde¬ 
pendent works. Ref. 30 pointed out that the FQH states 
described by wave functions 

= [Xn{{zi})]'^, 

or ^u=n/(m+n)(Xzi}) = Xl{{Z‘i})[Xn{{Zi})] (2) 

have excitations with non-abelian statistics, where Xn 
is the fermion wave function of n-filled Landau lev¬ 
els. The edge of the above FQH states are described 
by U{l)'^"^/SU{m)n or t7(l)"’'"+®/S'C/(m)„ Kac-Moody 
current algebra.®®^®® Those results were obtained by de¬ 
riving their low energy effective SU(m) level n Chern- 
Simons theory or 17(1) x SU{m) level n Chern-Simons 
theory. In the same year, Ref. 35 conjectured that the 
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FQH state described by p-wave paired wave function^®’^^ 





1 


LZi -Z2Z3- Z4 




(3) 


has excitations with non-abelian statistics. Its edge 
states were studied numerically in Ref. 38 and were found 
to be described by a c = 1 chiral-boson conformal field 
theory (CFT) plus a c = 1/2 Majorana fermion CFT. 
Such a result about the edge states supports the conjec¬ 
ture that the p-wave paired FQH state is non-abelian, 
since the edge for abelian FQH states always have inte¬ 
ger chiral central charge A few years later, the 

non-abelian statistics in p-wave paired wave function was 
also confirmed by its low energy effective SO{5) Chern- 
Simons theory.^® 

It is interesting to point out that long before the dis¬ 
covery of FQH states, Onnes discovered superconduc¬ 
tor in 1911.^^ The Ginzburg-Landau theory for symme¬ 
try breaking phases is largely developed to explain su¬ 
perconductivity. However, the superconducting order, 
that motivates the Ginzburg-Landau theory for symme¬ 
try breaking, itself is not a symmetry breaking order. 
Superconducting order (in real life with dynamical C/(l) 
gauge field) is an order that is beyond Landau symme¬ 
try breaking theory. Superconducting order (in real life) 
is an topological order (or more precisely a Z 2 topologi¬ 
cal order).jg quite amazing that the experimental 
discovery of superconducting order did not lead to a the¬ 
ory of topological order, but instead, lead to a theory of 
symmetry breaking order, that fails to describe super¬ 
conducting order itself. 


II. WHAT IS TOPOLOGICAL ORDER? 

A. Topological ground state degeneracy 

The above description of topological order is highly 
incomplete and highly unsatisfactory. This is because the 
characterization of topological order is through specifying 
what it is not: topological order is a kind of orders that 
cannot be described by symmetry breaking. But what is 
the topological order? 

To appreciate the difficulty of describing topological 
order, let me tell a story about a tribe. The tribe uses a 
language that contains only four words for counting: one, 
two, three, and many-many. It is very hard for a tribe 
member to describe a naturally occurring phenomenon ~ 
a large herd of deers. He can only describe the number 
of deers in the herd by what it is not - the number is not 
one, nor two, nor three. 

Similarly, the possible organizations of many particles 
in naturally occurring states can be very rich, much richer 
than those described by symmetry breaking. To describe 
the new orders (such as the topological orders), we need 
to introduce new tools and new languages. The richness 


of nature is not bounded by the known theoretical for¬ 
malism. The Landau’s symmetry breaking theory corre¬ 
sponds to “one”, “two”, “three” which describes a small 
class of orders. Many other orders also exist in nature, 
but we do not know how to describe them. Therefore, we 
introduced terms like “spin liquid”, “non-Fermi liquid”, 
“exotic order”, “preformed pair”, “dynamical stripe”, 
etc . Just like the term “many many” in the above story, 
those terms mainly describe what it is not than what it 
is. 

The symmetry breaking theory is the only language 
that we know to describe phases and orders. But topo¬ 
logical order, by definition, cannot be described by the 
symmetry breaking theory. If we abandon the only lan¬ 
guage that we know, how can we say any thing? Where 
do we start to understand the topological order? So the 
development of topological order theory is mainly trying 
to come up with a proper way to name/label topologi¬ 
cal orders. We hope the name/label to carry information 
that allows us to derive all the universal properties of the 
corresponding topological order from its name/label. 

To make progress, let us point out that, in physics, to 
define and to introduce a concept is to design an exper¬ 
iment (a laboratory one or a numerical one). We need 
to identify measurable quantities such that the measure¬ 
ment of those quantities facilitate the definition of the 
concept. So in physics, once you design an experiment, 
you define a concept. And only after you design an ex¬ 
periment, do you define a concept. 

So what experiments or what measurable quantities 
define the concept of topological order? It was noted that 
a. iz = Ijm Laughlin FQH state has m fold degenerate 
ground states on torus and a non-degenerate ground state 
on sphere.However, the different degeneracies was 
regarded as finite size and/or group theoretical effects 
without thermodynamical implications. 

In Ref. 22,23,29, it was shown that the ground state 
degeneracy of a chiral spin state or a FQH state is stable 
against any local perturbations, including random per¬ 
turbations that break all the symmetries.Thus the 
topology-dependent ground state degeneracies are a ro¬ 
bust or universal property with important thermodynam¬ 
ical implications: the topology-dependent and topologi¬ 
cally robust degeneracies can be used to define a phase (or 
a universality class) of a thermodynamical system {i.e. a 
system with a large size). So the topology-dependent 
ground state degeneracies is just what we are looking 
for: the measurable quantities (in a numerical experi¬ 
ment) that can be used to (partially) define topological 
order in chiral spin states and FQH states.Such kind 
of universal properties are also call topological invariants, 
since they are robust against any local perturbations. 

We would like to remark that the ground state degen¬ 
eracy discussed above is only an approximate degeneracy 
for a finite system, i.e. there is a small energy splitting 
e between different degenerate ground states. The en¬ 
ergy gap to other excited states is given by A (see Fig. 
I). It was shown in Ref. 23 and Ref. 29 that, for chi- 
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\ _ Nearly degenerate 

6 { = ground states 

FIG. 1: The energy levels of a topologically ordered state 
with a finite size. The splitting e of the nearly degenerate 
ground states approaches to zero in the large-system limit. 


ral spin states and FQH states, e is exponentially small: 
e ^ while A is finite in the limit where the system 

size approaches infinite: L —> oo. 

The topology-dependent ground state degeneracy is 
an amazing phenomenon. In both FQH and chiral spin 
states, the correlation of any local operators are short 
ranged. This seems to imply that FQH and chiral spin 
states are “short sighted” and they cannot know the 
topology of space which is a global and long-distance 
property. However, the fact that ground state degen¬ 
eracy does depend on the topology of space implies that 
FQH and chiral spin states are not “short sighted” and 
they do find a way to know the global and long-distance 
structure of space. So, despite the short-ranged correla¬ 
tions of all the local operators, the FQH and chiral spin 
states must contain certain hidden long-range structure. 
The robustness of the ground state degeneracy suggests 
that the hidden long-range structure in FQH/chiral-spin 
states is also robust and universal. A term topological 
order was introduced to describe such a “robust hidden 
long range structure” 

More recently, such a “robust hidden long range 
structure” was identified to be the long-range entan¬ 
glement defined by local unitary transformations.®^^®"* 
Thus topological order is nothing but the pattern of long 
range entanglement. Different patterns of long-range en¬ 
tanglement (or different topological orders) correspond 
to different quantum phases. Chiral spin liquids, 
integral/fractional quantum Hall states®’^®’®®, Z 2 spin 
liquids, ®®“®® non-Abelian FQH states,®®’®®’®®’®® etc are 
examples of topologically ordered or long-range entan¬ 
gled phases. 


B. Topological order and phase transitions 

In section I A, we define a quantum phase as a region 
bounded by lines of singularity in the ground state en¬ 
ergy (or some other local quantities). In section H, we 
define a topologically ordered phase as a region character¬ 
ized by a certain ground state degeneracy. Are these two 
definition self consistent? As one topologically ordered 
phase changes into another topologically ordered phase, 
the ground state degeneracy may change from one value 



(a) 



FIG. 2: Two types of phase transitions between two gapped 
states, (a) Energy levels of a Hamiltonian Ha{g) as functions 
of the coupling constant g. The topologically order state for 
g < gc changes into a trivial state for g > gc via a first order 
phase transition, (b) Energy levels of a Hamiltonian Hb{g) as 
functions of the coupling constant g. The topologically order 
state for g < gc changes into a trivial state for g > gc via a 
continuous phase transition. In this case, the ground state of 
Hb{gc) is a quantum critical state. 



FIG. 3: An impossible g dependence of energy levels. 


to another value. So why a change in the ground state 
degeneracy corresponds to a singularity in some averages 
of local quantities? 

We note that the ground state degeneracy that charac¬ 
terize topological order is robust again any perturbations. 
So a small change in the Hamiltonian will not change 
the ground state degeneracy. However, a large change of 
Hamiltonian can cause a change in ground state degen¬ 
eracy. The ground state degeneracy can change in two 
different ways as described by Fig. 2. 

In Fig. 2a, the ground state degeneracy changes due 
to an level crossing. The ground state energy has a dis¬ 
continuous first order derivative at the crossing point. 
The corresponding phase transition is a first order phase 
transition. 

Since the ground state degeneracy is robust against any 
perturbations, this means that the degenerate ground 
state wave functions are locally indistinguishable. The 
ground states cannot be split without losing the locally 
indistinguishable property. The only way to lose locally 
indistinguishable property is to develop a long range cor¬ 
relation i.e. closing the energy gap of other excitations. 
So, the situation described by Fig. 3 cannot happen. The 
possible situation is described by Fig. 2b, where the en- 
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ergy gap of the excitations closes as g —)■ and reopens 
as g passes Qc- The closing of the energy gap allow the 
ground state degeneracy to change. The closing of the 
energy gap at gc cause a singularity in the ground state 
energy (or some other local quantities). Such a phase 
transition is a continuous phase transition. 

We see that the change of ground state degeneracy 
of topologically ordered state and singularity in ground 
state energy always happen at the same place. Thus the 
topological ground state degeneracy characterize a phase 
and a change of the topological ground state degeneracy 
marks a phase transition. 


C. Topological invariants — Towards a complete 
characterization of topological orders 

Soon after the introduction of topological order 
through topologically robust and topology-dependent 
ground state degeneracies, it was realized that the 
topology-dependent degeneracies are not enough to char¬ 
acterize all different topological orders [note that, as dis¬ 
cussed in section I A, orders are defined through phase 
transitions]. Certain different topological orders can have 
exactly the same set of ground state degeneracies for all 
compact spaces. 

To obtain a complete topological invariant that can 
fully characterize topological orders, in Ref. 23,61, it 
was conjectured that the non-Abelian geometric phases^^ 
(both the U{1) part and the non-Abelian part) of degener¬ 
ate ground states generated by the automorphism of Rie- 
mann surfaces can completely characterize different topo¬ 
logical order sf^ 

We note that an automorphism of a Riemann surface 
change the Hamiltonian H defined on the surface to an¬ 
other Hamiltonian H' which is defined on the same sur¬ 
face. If we smoothly deform the Hamiltonian H to H', 
plus the automorphism transformation at the end, we 
will get a family of Hamiltonians that form a “closed 
loop” ,23.61 yyg ^gg gggj^ g loop-like deformation path 
of Hamiltonians, with their degenerate ground states, to 
define a non-Abelian geometric phase.Thus, for every 
automorphism of Riemann surfaces, we can produce a 
non-Abelian geometric phase which is a unitary matrix. 

Such a unitary matrix is uniquely determined by the 
automorphism (up to a path dependent over all 17(1) 
phase). To understand such a result, let us assume that 
the unitary matrix is not uniquely determined by the 
automorphism, i.e. a small change of deformation path 
leads to a different unitary matrix beyond the different 
over all 1/(1) phase. This will mean that the small change 
of deformation path causes different phase shifts for dif¬ 
ferent degenerate ground states. Since the small change 
of deformation path are local perturbations, the differ¬ 
ent phase shifts for different degenerate ground states 
will implies that the degenerate ground states are locally 
distinguishable, which contradict the robustness of the 
degeneracy against any local perturbations and the lo¬ 


cally indistinguishable property of the degenerate ground 
states. 

As a result, the above unitary matrices form a projec¬ 
tive representation of automorphism group of Riemann 
surfaces.The automorphism group GAut contain 
a connected subgroup G^ut- GAut/G^ut is the mapping 
class group (MCG). We note that the non-Abelian geo¬ 
metric phases for the automorphisms in G^^^ are all pure 
17(1) phases, since the loops that correspond to the auto¬ 
morphisms in G5 ,^^j are all contractible to a trivial point. 
Thus the non-Abelian geometric phase also generate a 
projective representation of MCG. 

We see that the non-Abelian geometric phases contain 
a universal non-Abelian part^^’®^ and a path dependent 
Abelian part^^’®^. The non-Abelian part carries infor¬ 
mation about the projective representation of MCG. For 
torus, the MCG is SL{2,Z), which is generate by a 90° 
rotation and a Dehn twist. For such two generators of 
MCG, the associated non-Abelian geometric phases is de¬ 
noted by S and T, which are unitary matrices. S and T 
generate a projective representation of MCG SL{2, Z) for 
torus. 

The Abelian part of the non-Abelian geometric 
phases is also important: it is related to the gravita¬ 
tional Chern-Simons term®^’®®^®® and carries informa¬ 
tion about the chiral central charge c for the gapless 
edge excitations.®^’®® The chiral central charge c can 
be measured directly via the thermal Hall conductivity 

Kr = cf of the sample.®®’®^ 

It is believed that (S', T, c) form a complete and one- 
to-one description of 2-l-lD topological orders, which is 
consistent with the previous conjecture in Ref. 23. So 
(S, T, c) are the new words, like “four”, “five”, “six”, 
“ten”, “eleven”, “twelve”, etc in our tribe story, that 
we are looking for, to describe/label topological orders. 
Since {S,T,c) may completely describe 2-l-lD topolog¬ 
ical orders, we may be able develop a theory of 2+lD 
topological order based (S, T, c). 


D. Wave-function-overlap approach to obtain S,T 

We like to mention that in addition to use non-Abelian 
geometric phases to obtain (S, T) matrices, there are sev¬ 
eral other ways to obtain them.®®^^^ In particular, one 
can use wave function overlap to extract (S, T) matrices 
directly from the degenerate ground states wave func¬ 
tions on torus, provided that the system have translation 
symmetry.®^’^®^^® (The non-Abelian-geometric-phase ap¬ 
proach can obtain [S, T) matrices even from systems 
without translation symmetry.) It was argued that for a 
system on a d-dimensional torus of volume V with the 

set of topologically degenerate ground states 
the overlaps of the degenerate ground states have the 
following form 


( 4 ) 
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where W are transformations of the wave functions in¬ 
duced by the MCG transformations of the space 
T^, / is a non-universal constant, and is an univer¬ 
sal unitary matrix. 

We know that a MCG transformation W maps the 
space T'^ to itself: —>■ It transforms a ground state 

wave function on space to another wave function 
W|'0j) on the same space T^. Since the MCG transfor¬ 
mation W is not a symmetry of the Hamiltonian, the 
new wave function W\ipj) is not longer a ground state of 
the Hamiltonian. So the overlap of W\il}j) with a ground 
state \'ipi) is exponentially small in large volume limit. 
('0i|IH|'0j) ~ It seems that such an overlap con¬ 

tains no useful universal information about topological 
order. What was discovered in Ref. 64 is that if we sep¬ 
arate out the volume dependent exponential factor, the 
volume-independent constant factor contains useful 
universal information about topological order. 

We note that the volume-independent constant factor 
is a unitary matrix. In contrast to non-Abelian ge¬ 
ometric phases, such a unitary matrix has no C/(I) phase 
ambiguity. Those unitary matrices (from different MCG 
transformations W) form a representation of the MCG 
of the space MCG(T‘^) = SL{d,Z), which is robust 
against any perturbations. For 2-l-lD cases, the MCG of 
the torus is generate by 90° rotation S and Dehn twist 
T. The corresponding unitary matrices S —>■ = S 

and T —>■ M'^ = T generate a unitary representation of 
SL{2, Z) (instead of a projective representation as for the 
case of non-Abelian geometric phases). As a result, we 
can use a unitary representation of MCG, {S,T), plus 
the chiral central charge c to characterize all the 2-l-lD 
topological orders. 

We also like to point out that we can always choose a 
so called excitation basis for the degenerate ground state 
(see Section D). In such a basis, T is diagonal and Su 
are real and positive. It is {S, T) in such a basis, plus the 
chiral central charge c, that may fully characterize all the 
2-1-lD topological orders. 


E. The current systematic theories of topological 
orders 


integer AT-matrices to classify all of them.^® So the inte¬ 
ger matrices K are also the new words, like “4”, “5”, “6”, 
etc in our tribe story, that can be used to describe/label a 
subset of topological orders - Abelian topological orders. 
Such a AT-label completely determine the low energy uni¬ 
versal properties of the corresponding topological order. 
For example, the low energy effective theory for the topo¬ 
logical order labeled by K is given by the following ?7(1) 
Chern-Simons theory^^’^®’®^^®® 

C = ^ai^d,ajxe>^’'\ (5) 

47r 

Such an effective theory or the topological order labeled 
by K can be realized by a concrete physical system - a 
multi-layer FQH state: 


n ~ 4n ~ ® 

I]i<j 




( 6 ) 


where z( = xj + ij/f is the coordinate of the particle 
in /**' layer. 

Gertainly, the topological order described by K are also 
described by {S, T, c). The non-Abelian geometric phases 
for some canonical choices of path are calculated for the 
bosonic K topological order described by eqn. (6):®^ 


(7) 

where c is the difference in the numbers of positive and 
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negative eigenvalues of K, 
two times the spin vector of the Abelian FQH state, and 
= (iVi,7V2,... ,1V„) are the numbers of the bosonic 
“electrons” in each layer.®® We see that the 17(1) factors 
depend on the number of electrons and are not univer¬ 
sal. But we can isolate the universal non-Abelian part 
by taking the limit AT —>■ 0, and find 


T i i2777^ ^ir 

«/3 = e 24 e dot(3 

g— i2-Kf3^ K(x 


\/|det(77)| 


( 8 ) 


We see that the universal non-Abelian part of the non- 
Abelian geometric phases determines S', T, and c mod 


24. 


We like to remark that topological order {i.e. long- 
range entanglement) is truly a new phenomena. They 
require new mathematical language to describe them. 
Some early researches suggest that tensor category 
theory®®7®“®® and simple current algebra®^’®®’®^’®® (or 
pattern of zeros®®^®®) may be part of the new mathe¬ 
matical language. Using tensor category theory, we have 
developed a systematic and quantitative theory that clas¬ 
sify topological orders with gappable edge for 2-l-lD in¬ 
teracting boson and fermion systems. 78.so 

For 2-1-ID topological orders (with gapped or gapless 
edge) that have only Abelian statistics, we have a more 
complete and simpler result: we find that we can use 


III. TOPOLOGICAL EXCITATIONS 

We have seen that we can use unitary representation of 
MGG and the chiral central charge, (S', T, c), to charac¬ 
terize/label/name all the 2-1-ID topological orders. It is 
possible that (S, T, c) is a full characterization of 2-l-lD 
topological orders, in the sense that all other universal 
properties of topological orders can be determined from 
the data (S, T, c). In this section, we will discuss some 
other universal properties of 2-l-lD topological orders, 
and see how those universal properties are determined 
by the data (S, T, c). 







excitation 

energy density\ ^ ground state 



FIG. 4: The energy density distribution of a quasiparticle. 

A. Local excitations and topological excitations 

Topologically ordered states in 2+ID are characterized 
by their unusual particle-like excitations which may carry 
fractional/non-Abelian statistics. To understand and to 
classify particle-like excitations in topologically ordered 
states, it is important to understand the notions of local 
quasiparticle excitations and topological quasiparticle ex¬ 
citations. 

First we define the notion of “particle-like” excitations. 
Consider a gapped system with translation symmetry. 
The ground state has a uniform energy density. If we 
have a state with an excitation, we can measure the en¬ 
ergy distribution of the state over the space. If for some 
local area, the energy density is higher than ground state, 
while for the rest area the energy density is the same as 
ground state, one may say there is a “particle-like” exci¬ 
tation, or a quasiparticle, in this area (see Figure 4). 

Quasiparticles defined like this can be divided into two 
types. The first type can be created or annihilated by 
local operators, such as a spin flip. So, the first type of 
particle-like excitation is called local quasiparticle excita¬ 
tions. The second type cannot be created or annihilated 
by any finite number of local operators (in the infinite 
system size limit). In other words, the higher local en¬ 
ergy density cannot be created or removed by any local 
operators in that area. The second type of particle-like 
excitation is called topological quasiparticle excitations. 

From the notions of local quasiparticles and topological 
quasiparticles, we can further introduce the notion topo¬ 
logical quasiparticle type, or simply, quasiparticle type. 
We say that local quasiparticles are of the trivial type, 
while topological quasiparticles are of nontrivial types. 
Two topological quasiparticles are of the same type if 
and only if they differ by local quasiparticles. In other 
words, we can turn one topological quasiparticle into the 
other one of the same type by applying some local oper¬ 
ators. 


B. Fusion space and internal degrees of freedom 
for the quasiparticles 

The quasiparticles have locational degrees of freedom, 
as well as internal degrees of freedom. 

To understand the notion of internal degrees of free¬ 
dom, let us discuss another way to define quasiparticles: 
Consider a gapped local Hamiltonian qubit system de¬ 
fined by a local Hamiltonian Hq in d dimensional space 
M‘^ without boundary. A collection of quasiparticle ex¬ 
citations labeled by i and located at Xi can be produced 


as gapped ground states of Hq -\- 5H where SH is non¬ 
zero only near Xi ’s. By choosing different 6H we can 
create (or trap) all kinds of quasiparticles. We will 
use ii to label the type of the quasiparticle at Xi. The 
gapped ground states of Hq -|- 5H may have a degen¬ 
eracy II(M^; ii, * 2 ) • • •) which depends on the quasipar¬ 
ticle types ii and the topology of the space M^. The 
degeneracy is not exact, but becomes exact in the large 
space and large particle separation limit. We will use 
V{M‘^;ii,i 2 ,- ‘' ) to denote the space of the degenerate 
ground states. If the Hamiltonian Hq-\-5H is not gapped, 
we will say * 2 , • ’ ’) = 0 (i-e-, H, * 2 , • • •) 

has zero dimension). If Hq -\- 6H is gapped, but if SH 
also creates quasiparticles away from xfs (indicated by 
the bump in the energy density away from xfs), we will 
also say D{M‘^; * 2 , •' ’) = 0- (In this case quasiparti¬ 

cles at xfs do not fuse to trivial quasiparticles.) So, if 
H, 12 , • • •) >0, SH only creates/traps quasiparti¬ 
cles at XiS. 

If we choose the space to be a d-dimensional sphere 
Md — then the number of the degenerate ground 
states, D{S'^;ii,i 2 , ■'') represents the total number 
of internal degrees of freedom for the quasiparticles 
fi,r 2 ,---). To obtain the number of internal degrees 
of freedom for type-i quasiparticle, we consider the di¬ 
mension D{S^] i,i, - ■ ■ ,i) of the fusion space on n type-i 
particles on S'^. In large n limit D{S‘^\i,i,- ■ ■ ,i) has a 
form 

lnD{S'^;i,i, ■ ■ ■ A) = n{\Tidi -I- o(I/n)). (9) 

Here di is called the quantum dimension of the type- 
i particle, which describe the internal degrees of free¬ 
dom the particle. For example, a spin-0 particle has a 
quantum dimension d = 1, while a spin-1 particle has 
a quantum dimension d = 3. For particles with abelian 
statistics, their quantum dimensions are always equal to 
1. For particles with non-abelian statistics, the quantum 
dimensions d > 1, but in general the quantum dimensions 
d may not be integers. 


C. Simple type and composite type 

Even after quotient out the local quasiparticle excita¬ 
tions, topological quasiparticle type still have two kinds: 
simple type and composite type. 

We can also use the traping Hamiltonian Hq -\- SH and 
the associated fusion space V(M‘^; * 2 , • • •) to under¬ 

stand the notion of simple type and composite type. 

If the degeneracy * 2 , ■ ’ ’) (the dimension of 

* 2 , • • •)) cannot not be lifted by any small local 
perturbation near Xi, then the particle type H at Xi is 
said to be simple. Otherwise, the particle type ii at Xi 
is said to be composite. 

When ii is composite, the space of the degenerate 
ground states V(M‘^; ii, ^ 2 , is, • • •) has a direct sum de- 
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composition: 

= VCM''; ji, i2, ^3, • • •) ® V(M® fci, Z2, is, • • •) 

0 V(M‘^;/i, i2, Z3, • • •) 0 • ■ ■ (10) 

where ji, /ci, Zi, etc. are simple types. To see the 
above result, we note that when ii is composite the 
ground state degeneracy can be split by adding some 
small perturbations near Xi. After splitting, the original 
degenerate ground states become groups of degenerate 
states, each group of degenerate states span the space 
or fci, Z 2 , 13 , • • •) etc. which 

correspond to simple quasiparticle types at Xi. The 
above decomposition allows us to denote the composite 
type as 

*1 = ji 0 fci 0 Zi 0 ■ ■ • . (11) 


- ■ ■) is determined completely by the fusion 

rules Nl^. 

Let us then consider the fusion of 3 simple quasipar¬ 
ticles i,j,k. We may first fuse i,j, and then with fc, 
{i(^j)^k = 0 fc = 0i(E„^ lVbA^™fc)Z. We 

may also first fuse j, k and then with i, i ® (j ® k) = 
i ® = 0i(Em^r^m )^- The two ways of 

fusion should produce the same result and this requires 
that 

= (14) 

m m 

Note that here, we do not require . 

The fusion coefficients are also topological invari¬ 

ants of the topological order. (S', T, c) can determine such 
topological invariants. In fact, S alone can determine 


The degeneracy Z1(M‘^; ii, Z 2 , • • •) for simple particle 
types ii is a universal property (i.e., a topological invari¬ 
ant) of the topologically ordered state. In this paper, 
when we said particle/topological type, we usually mean 
simple type. The number of simple types (including the 
trivial type) is also a topological invariant of the topo¬ 
logical order. Such a number is referred as the rank of 
the topological order. 

We have claimed that (S, T, c) can determine all other 
topological invariants of a topological order, including its 
rank. Indeed, the dimension of the S or T matrices is the 
rank of the topological order. 


< = 1: 


SuSi,{SikY 

Sn 


(15) 


which is the famous Verlinde formula.®^ 

The internal degrees of freedom {i.e. the quantum di¬ 
mension di) for the type-i simple particle can be calcu¬ 
lated directly from . In fact di is the largest eigen¬ 
value of the matrix Ni, whose elements are {NY^j = . 

We see that S matrix determines the internal degrees of 
freedom of the simple particles. 


D. Fusion of quasiparticles 

When we fuse two simple types of topological particles 
i and j together, it may become a topological particle of 
a composite type: 

i 0 j = Z = fci 0 ^2 0 • • • , (12) 

where i,j, ki are simple types and Z is a composite type. 
Here, we will use an integer tensor to describe the 
quasiparticle fusion, where i,jjk label simple types. Such 
an integer tensor is referred as the fusion coefficients 
of the topological order, which is a universal property of 
the topologically ordered state. 

When = 0, the fusion of i and j does not contain 
k. When = I, the fusion of i and j contain one k: 
i05 = fc0Zci0Ze2 0- -- . When = 2, the fusion of i 
and j contain two fc’s: z 0 j = fc 0 Zc 0 fci 0 ^2 0 • • • ■ This 
way, we can denote that fusion of simple types as 

i®j = ®kNYk. ( 13 ) 

In physics, the quasiparticle types always refer to sim¬ 
ple types. The fusion rules is a universal prop¬ 
erty of the topologically ordered state. The degeneracy 


E. Quasiparticle intrinsic spin 

For 20lD topological orders, the quasiparticles can 
also braid. We also need data to describe the braiding of 
the quasiparticles in addition to the fusion rules We will 
discuss the braiding in this and next subsections. 

If we twist the quasiparticle at Xi by rotating 5H at a;i 
by 360° (note that 5H aixi has no rotational symmetry), 
all the degenerate ground states in V(M‘^; * 2 , * 3 , • • •) 

will acquire the same geometric phase e'^u provided that 
the quasiparticle type ii is a simple type. This is because 
when ii is a simple type, no local perturbations near Xi 
can split the degeneracy. Thus the degenerate ground 
states are locally indistinguishable near Xi. As a result, 
the 360° rotation cause the same phase shift e‘®u for all 
the degenerate ground states. We will call Si = |^ mod 
I the intrinsic spin (or simply spin) of the simple type z, 
which is another universal property of the topologically 
ordered state. (5, T, c) can determine the topological in¬ 
variants Si as well. In fact, Si mod I are given by the 
eigenvalues or the diagonal elements of T and c: 

( 16 ) 

(note that T is diagonal in the excitation basis). 
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F. Quasiparticle mutual statistics 


If we move the quasiparticle *2 at around the quasi¬ 
particle zi at Xi, we will generate a non-Abelian geo¬ 
metric phase - a unitary transformation acting on the 
degenerate ground states in V{M‘^; is, - ■ ■). Such a 
unitary transformation not only depends on the types ii 
and Z 2 , but also depends on the quasiparticles at other 
places. So, here we will consider three quasiparticles of 
simple types i, j, fc on a 2D sphere . The ground 
state degenerate space is V{S^-,i,j,k). For some choices 
of i, j, k, j,k) > 1, which is the dimension of 

j,k). Now, we move the quasiparticle j around 
the quasiparticle i. All the degenerate ground states in 
k) will acquire the same geometric phase 


e 


qW 


giSTTSfc 

0i27rsi Qi2nsj 


(17) 


This is because, in j,k), the quasiparticles i 

and j fuse into k (the anti-quasiparticle of k). Moving 
quasiparticle j around the quasiparticle i plus rotating 
i and j respectively by 360° is like rotating k by 360°, 

i.e. e'^'J Qi^Tvsi ^i 2 Trsj _ gi27rsfc^ This leads to eqn. (17). 
We see that the quasiparticle mutual statistics is deter¬ 
mined by the quasiparticle spin Si and the quasiparticle 
fusion rules . For this reason, we call the set of data 
{Nl\si) quasiparticle statistics. 

In fact, in order for data to describe a valid 

quasiparticle statistics, they must satisfy certain condi¬ 
tions (known as Vafa theorem®®“^°°). Let us consider 
the fusion space V{S‘^]i,j,k,l). Let Wij be the non- 
abelian geometric phase {i.e. the unitary matrix acting 
V{S^-,i,j,k,l)) generated by moving particle i around 
particle j, Wi^k by moving particle i around particle k, 
and Wijk by moving particle i around both particle j 
and k (see Fig 5). We see that Wijk = tFqfetFij, or 


det{Wijk) = dei{Wi^k) det{Wij) 


We note that 


det(IFij) = n ( 

r 

det(Wi,/c) = n ( 

r 

det(Wijfe) = n ( 


gi27rsr 

^ 1 2 TT S ^ ^ 1 2 TT S j 


Aff Afp 


gi27rsf . ATpATp 

gi27rsi 0i27rs.p J 


. 1 27rSr 


) 


This way, we obtain 

g i 27r Sr- (Afp Arp -I-Arp ATp -I-ATp ATp) 

_ „i27ry:^ Si(7VpArp-ArpArf'=-ArpA/p) ^ 

— 0 ^ t I i X 

Arp) 

_ i27r(s^+Sj -|"Sfc+Si) TVpTVP 


(18) 


(19) 


(20) 


where the properties eqn. (29) and eqn. (38) are used. 


Wijk 



FIG. 5: The braiding procedure to derive Vafa theorem. 


IV. A THEORY OF 2+lD BOSONIC 
TOPOLOGICAL ORDERS 


A. A theory of 2+lD topological orders based on 

(s,r,c) 

In this section, we would like to develop a theory of 
2-I-ID topological orders based on {S, T, c). We have seen 
that we can measure {S,T,c) for every 2-l-ID topologi¬ 
cal orders (in particular using the wave function over¬ 
lap eqn. (4)), and every 2-l-lD topological orders are de¬ 
scribed by {S, T, c) where S, T are unitary matrices and 
c is a rational number. However, not every {S, T, c) can 
describe existing topological orders in 2-l-lD. So to de¬ 
velop a theory of topological order based on (S, T, c), we 
need to find the conditions on {S, T, c). If we find enough 
conditions on {S,T,c), then every {S,T,c) that satisfies 
those conditions will describe an existing topological or¬ 
der. This way, we will have a theory of topological orders. 

So here, we will follow Ref. I01-I04 and list the known 
conditions satisfied by a {S, T, c) that corresponds to an 
existing 2-l-lD topological order: 

{S, T, c) conditions: 

1. S is symmetric and unitary with Pn > 0, and sat¬ 
isfies the Verlinde formula:®^ 

ivp= (21) 

where i,j,--- = 1,2, ■ ■ ■ , N. IVp is called fusion 
coefficient, which gives the fusion rule for quasi¬ 
particles. 

2. Let 

d^ = ^ (22) 

*11 

which is called quantum dimension. Then di > 1 
is the largest eigenvalue of the matrix N^, whose 
elements are {Ni)kj = . 

3. T is unitary and diagonal: 

Ty = (23) 

Here Si is called topological spin {9i = 27rsi is called 
statistical angle), c is the chiral central charge. 







11 


4. S and T satisfy: 

{STf = S^ = C, C^ = l, Cij=N\^. (24) 

Thus S and T generate a unitary representation of 
SL{2,Z). 

5. S and T also satisfy [see eqn. (223) in Ref. 105]: 

(25) 

k 

where D = i/X)* d'f. 

6. N]^ and also satisfy®®^^™ (see eqn. (20)) 


l.fV^ are non-negative integers that satisfy 

k^l 

n n 

Y, = Y NrNi'^ or NkN, = N,Nk (29) 

m—1 n—1 

where = 1, 2, • • • , n, and the matrix Ni is 

given by {Ni)kj = ■ In fact Nl^ defines a charge 

conjugation i ^ i: 

Nl^ = 6-^. (30) 


Y ^^kiSr = 0 mod 1 


(26) 


We also refer n as the rank of the corresponding 
topological order. 


where 

V[,ki = 


- {5^r + 5jr + 5kr + Sir) Y 

ra 


(27) 


7. Let 

(28) 

jk 

then^°®d04 = Q jf j ^ jy^ = j-j if J = 1. 

The above are the necessary conditions in order for 
(S', T, c) to describe an existing 2-l-lD topological order. 
In other words, the (S, T, c)’s for all the topological or¬ 
ders are included in the solutions. 

However, it is not clear if those conditions are sufh- 
cient. So it is possible that some solutions are “fake” 
(S, T, c) that do not correspond to any valid topological 
order. It is also possible that a valid solution (S, T, c) 
may correspond to several topological orders. 

To see if there are any “fake” (5”, T, c)’s in our lists. 
Ref. 106 tries to construct explicit many-body wave func¬ 
tions for those (S', T, c)’s in the lists, using simple current 
algebra.®"'’®^’®^ We find that all the (S, T, c)’s in our lists 
are valid and correspond to existing topological orders. 


2. and Si satisfy®® (see eqn. (20)) 

Y ^^jkiSr = 0 mod 1 (31) 


where 

- {S„ + <5,, + Skr + Sir) Y (32) 

m 

Those are the conditions that allows us to show Si 
and c to be rational numbers. 


3. Let di be the largest eigenvalue of the matrix W- 
Let 


S^J = 


1 


—(33) 

V 2^i k 
Then, S is unitary and satisfies®^ 

Sii>0, Nl^ = ^ ^i-^hi^ik)* 


4. Let 


B. 


A theory of 2+lD topological orders based on 


{K',s.,c) 


T = e 


277 i Si ^ — 277 i 7^ 


24 A - • 

^7,7 • 


(35) 


Then 


From the above conditions, we see that, instead of us¬ 
ing (5,r,c), we can also use (iV^®, Sj, c) to describe topo¬ 
logical orders, since (S', T, c) can be expressed in terms of 
, Si, c), and ,Si,c) can be expressed in terms of 
{S,T,c). So we can develop a theory of topological or¬ 
ders based on ,Si, c), instead of (S, T, c). Again not 
all (Nl^Si, c) describe existing 2-l-lD topological orders. 
Here we list the necessary conditions on [N]^ ,Si,c): 

,Si,c) conditions: 


(ST)® = S^ = C, C® = 1. (36) 

In fact = N\K 

5. Let 

jk 

Then®®®d04 jy^ _ q ^ ^ ^ 






12 


The above conditions are necessary for {Nl^ ,Si, c) to de¬ 
scribe an existing 2-l-lD topological order. If the above 
conditions are also sufficient, then the above will repre¬ 
sent a classifying theory of 2-l-lD topological orders. 

In section V, we will solve the above conditions to ob¬ 
tain a list 2-1-ID topological orders. We like to men¬ 
tion that solving the above conditions is closely related 
to classifying modular tensor categories. Ref. 103 have 
classified all the 70 modular tensor categories with rank 
N = 1,2,3,4, using Galois group. In this paper, we will 
try to solve the above conditions numerically for higher 
ranks. 


V. 2+lD TOPOLOGICAL ORDERS WITH LOW 
RANKS AND LOW QUANTUM DIMENSIONS 

A. A numerical approach 


where D* is a diagonal matrix given by {D^)ii = d\ = 
Sii/Sii- We see that even though Ni may not be hermi- 
tian, we still require that 

Ni can be diagonalized by a unitary matrix (42) 

This is the third kind of conditions on 

To get more information, let ui be the common eigen¬ 
vectors of a set of W’s, i G I and / C {1, • • • , N}. We will 
try to calculate S from such a subset of N^s. Let U be a 
linear combination of the set of W’s, V = 

Let d] be the eigenvalue of W for the eigenvector ui. Let 
I belong to the set of indices that label eigenvectors that 
have non-degenerate eigenvalues for V. In this case, the 
corresponding eigenvector uj is unique up to a C/(I) phase 
factor. Then those non-degenerate normalized eigenvec¬ 
tors with the first element being positive satisfies 


Here, we will assume the conditions in section IV B 
to be sufficient, and treat them as a classifying theory of 
2-l-lD topological orders. In this section, we will describe 
how to numerically solve those conditions to obtain a 
list of simple 2-l-lD topological orders. Our approach is 
similar to that used in Ref. 102, where a list of fusion 
rings are obtain. Here, we will obtain a list of 2-l-lD 
bosonic topological orders. 

We first numerically solve the condition (1) in the 
{Nl^ ,Si,c) conditions in section IV B to obtain . Then 
we will use Smith normal form of integer matrices 

and/or Mij to solve the condition (2) to obtain a list of s^. 
We then use other conditions to obtain a list of , Si)’s 
that satisfy all those conditions by direct checking. The 
central charge c mod 8 is obtained from the condition 

(4)- 

To numerically solve the condition (1) in the , Si, c) 
conditions efficiently, it is important to find as many 
conditions on as possible. We first set Z = 1 in 
eqn. (29) and find the following symmetry condition on 

Nfj=Nj, = Nl = Nl. (38) 

The second kind of conditions on is that 

[N,,Nj]=0. (39) 


To find more conditions on , we note that since S 
is unitary, we may rewrite eqn. (34) as 




ShSi, 

Su 


or viN^ 


d]vi, 


(40) 


where the row eigenvector Vi is given by {vi)j = Sij and 
the eigenvalues dj = Su/Su. In other words there exist 
a symmetric unitary matrix that satisfies 


SNiS"^ = D\ (41) 


{u{),=SN^, (fo-)i^O. (43) 


where p is a permutation map I —>■ 1. For those (u[)j ’s, 
we have 




i ~ Q - 


(44) 


In other words ■ 

To summarize, let ui be the common eigenvectors of a 
set of Ni’s {i G I) with eigenvalue d], then 


(wf)i 7^ 0, 

3 


= 5. 


hh ’ 




* 


(45) 


for any i G I and /’s in the set of that label non¬ 
degenerate eigenvalues. If the above conditions are not 
satisfied, then corresponding Ni does not satisfy the nec¬ 
essary conditions to describe a topological order. 

Also, if the all the eigenvalues of V’s are non¬ 
degenerate, then ui determines S upto a permutation of 
the rows (see eqn. (43)). In this case, we can determine 
the full Nl^ using eqn. (21). 

We wrote a program to numerically search for ’s 
that satisfy the condition eqn. (38), eqn. (39), eqn. (45), 
and eqn. (21) (when all the eigenvalues of V are non¬ 
degenerate), by starting from {N^^ = Skj}, to {N^^, N^j}, 
to{N^^,N^^,N^^}, etc. 

After obtaining a list of fusion rules N]d , we then, for 
each fusion rule, use the Smith normal form of the integer 
matrix M to find sets of spins {si} that satisfy eqn. (??). 
Last, we select the combination {N^^ ,Si) that satisfy all 
the conditions and compute the central charge c in the 
process. This way we obtain a list of 2-l-lD topological 
orders. 
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TABLE I: A list of all 45 bosonic topological orders in 2+lD with rank = 1, 2, 3,4, 5 and with niax(A^^'’) < 3. The entries in 
red are the topological orders with niax(A^^'’) = 2. All other topological orders have = 0, 1. There are no topological orders 
with max(A^-’) = 3 and N < 5. The entries in blue are the composite topological orders that can be obtained by stacking lower 
rank topological orders. The first column is the rank N and the central charge c (mod 8). The second column is the topological 
entanglement entropy Stop = log 2 D, D — \/X]i Tli® quantum dimensions of topological excitations in the third column 
are expressed in terms of ■ The fourth column is the spins of the corresponding topological excitations. 



Stop 

di, d2,••• 

Si,S2,••• 

A^f 

‘S'top 

di, d2,• • • 

Sl, S2,• • • 

If 

0 

1 

0 





2f 

nB 

^14/5 

0.5 

0.9276 

1,1 

1,C3 

0,i 

o,i 

2fi 

2fl4/5 

0.5 

0.9276 

1,1 

bd 

0,-i 

0,-f 

3f 

3f/2 

oB 

'J3/2 

oB 

'J5/2 

oB 

'J7/2 

oB 

^8/7 

0.7924 

1 

1 

1 

1 

1.6082 

1,1,1 

1,1, d 

1,1, Cl 

1,1, Cl 

1,1, Cl 
i,C5Vcl 

0 i i 

3 ’ 3 

0 ^ ^ 

2’ 16 

0 i A 

2’ 16 

0 1 A 

2 ’ 16 

0 i ^ 

0-12 
’ 7 ? 7 

3f2 

oB 

• 2 - 1/2 

oB 

•2-3/2 

oB 

•2-5/2 

oB 

•2-7/2 

oB 

•2-8/7 

0.7924 

1 

1 

1 

1 

1.6082 

1 , 1,1 

1 , 1 , d 

1 , 1 , d 

1 , 1 , d 

1 , 1 , d 
i,d,cl 

0 -1 -1 

3 ’ 3 

2 ’ 16 

0 1 -A 

0 1 -A 

0 1 -A 

0 1 -2 
’ 7 ’ 7 

4 B.. 

4f 

4f 

4f9/5 

4® 

4-19/5 

4^ 

4i2/5 

4^ 

4i0/3 

1 

1 

1 

1 

1 

1.4276 

1.4276 

1.8552 

2.1328 

1,1,1,1 

1,1,1,1 

1,1,1,1 

1,1,1,1 

1,1,1,1 

1 , 1 , cl, cl 

1 , 1 , cl, cl 

1 , cl, cl, del 

i,d,d,C 7 

0 , 0 , 0 , 1 

0 , i,|,| 

0 1 1 1 

0 3 3 1 
^5 8 ’ 8 ’ 2 

0 1 1 1 

2 ’ 2 ’ 2 

0 i -A _2 

0 -1 A _2 

4 ’ 20 ’ 5 

0 -2 _2 1 

5 ’ 5 ’ 5 

0 i 2 _1 

3 ’ 9’ 3 

4-’" 

4fi 

4f2 

4f3 

41/5 

4® 

4i9/5 

4f’" 

4fl2/5 
4f 10/3 

1 

1 

1 

1 

1.4276 

1.4276 

1.8552 

1.8552 

2.1328 

1 , 1 , 1,1 

1 , 1 , 1,1 

1 , 1 , 1,1 

1 , 1 , 1,1 

1 , 1 , d,d 

1 , 1 , d, d 

1, d, d, dd 

1, d, d, dd 

l,d,C7fC7 

o,o,i,-i 

0-1-11 

0-3-3 1 

8 ’ 8 ’ 2 

0 -1 A 2 

4’ 20’ 5 

0 1 -A 2 
^’4’ 20’ 5 

0 1 —1 0 

5 ’ 5 ’ 

0 2 2 _1 

0-1-2 1 

^0 

rS.a. 

52 

rB,b 

^-2 

'^16/11 

^18/7 

1.1609 

1.7924 

1.7924 

2.5573 

2.5716 

1 , 1 , 1 , 1,1 

i,i,d,ci ,2 

i,i,d,ci ,2 

1 , cl, cl, cl, cl 

1, cl, cl, C 12 , C 12 

0 1 1 1 1 

’ 5 ’ 5 ’ 5 ’ 5 

n 2 2 1 5 

11 ’ 11 ’ 11 ’ 11 

0 1 113 

7 ? 7 ’ 7 ’ 7 

5f 

cB,b 

52 

rS,a 

^-2 

^-16/11 

^-18/7 

1.1609 

1.7924 

1.7924 

2.5573 

2.5716 

1,1,1,1,1 

i,i,Ci,d ,2 

i,i,Ci,d ,2 

bdiCdcdcl 

1 , cl, cl, C 12 , C 12 

0 2 2 2 2 

5 ’ 5 ’ 5 ’ 5 

0 0 -i ^ i 

8 ’ 8 ’ 3 

^2 2 15 

11 ’ 11 ’ 11 ’ 11 
Oil 1 3 

’ 7 ’ 7 ’ 7 ’ 7 


B. The stacking operation of topological order 


Before we present the result from the numerical cal¬ 
culation, let us discuss a stacking operation,®^ denoted 
by Kl. We note that stacking two rank N and rank 
N' topological orders described C = {Nl\si,c) and 
C = (A^'/ji ,s'^i,c') will give us a third topological order 
C” = Cm C with rank N" = NN' and 


4,=s, + 4, c" = c + c', 


d", = s" p = Stop + S^. 


(46) 


where Stop is the topological entanglement entropy 
Stop = logs 

The stacking operation Kl will make the set of topologi¬ 
cal order into a monoid. The trivial topological order Ctri 
(the product state) is the unit of the monoid. However, 
in general, a topological order C does not have an inverse 
respect to the stacking operation {i.e. there does not ex¬ 
ist a topological order C such that CMC = Cto). This 


is why the set of topological order only form a monoid 
instead of a group. However, some topological order does 
have an inverse respect to the stacking Kl operation. Such 
kind of topological orders are called invertible topological 

orders. ®®d09-ii2 

In 2-1-ID, the invertible topological orders form an 
Abelian group Z under to stacking Kl operation. The 
group is generated by the Ag FQH state described by 
the AT-matrix 




4100000 
1 2 1 0 0 0 0 
0 12 10 0 0 
0 0 12 10 0 
0 0 0 1 2 1 0 

0 0 0 0 1 2 1 

0 0 0 0 0 1 2 

Vo 0 0 0 1 0 0 


0 

0 

0 

1 

0 

0 

2/ 


(47) 


The Ag topological order C_Bs is invertible®® since it has 
no topological excitations®®’®^ (due to det{KEg) = 1). It 
is described by = (1,0,8). Stacking an Ag 



























TABLE II: A list of all 50 bosonic rank N — 6 topological orders in 2+ID with niax(A^^^) < 2. 




^top 




di , d 2 ,- 


Sl, S2,• • • 


comment 


6f 

6^1 

63® 

6^3 


1.2924 

1.2924 

1.2924 

1.2924 


1 , 1 , 1 , 1 , 1,1 


1 , 1 , 1 , 1 , 1,1 


_i 1 1 
12 ’ 12’ 4’ 3 ’ 3 

_i__ i_ 1 _i _i 

12 ’ 12 ’ 4’ 3’ 3 

111 

4 ’ 3 ’ 3 ’ 12 ’ 12 

_1 _1 _1 _5_ _5_ 

4’ 3 ’ 3 ’ 12 ’ 12 


2^1 Kl 3f 
2f Kl 3^2 


2f H 3 


2®i Kl 3^2 


O1/2 

rB 

0 - 1/2 

O3/2 

rB 

0-3/2 

O5/2 

0-5/2 

O7/2 

0-7/2 


1.5 

1.5 

1.5 

1.5 

1.5 

1.5 

1.5 

1.5 


1.1.1.1, d,C2 

1.1.1.1, Cl, Ci 
1,1,1,1, Cl, Cl 
1,1,1,1, Cl, Cl 
1,1,1,1, Cl, Cl 
1,1,1,1, Cl, Cl 
1,1,1,1, Cl, Cl 

1 , 1 , 1 , 1 , cl, cl 


1 

4 ’ ■ 
1 

4 ’ ■ 
1 

4 ’ ■ 
1 

4 ’ ■ 
1 

4 ’ ■ 
1 

4 ’ ■ 
1 

4 ’ ■ 
1 

4 ’ ' 


1 1 _3_ 

■ 4 ’ 2 ’ 16 ’ 16 

1 1 J_ 

■ 4 ’ 2 ’ 16’ 16 

1 i J_ A 

■ 4 ’ 2 ’ 16’ 16 
1 1 

■ 4 ’ 2 ’ 16 ’ 16 

1 1 _3_ X 

■ 4 ’ 2 ’ 16’ 16 
1 1 

■ 4 ’ 2 ’ 16 ’ 16 

1 i A _A 

■ 4 ’ 2 ’ 16’ 16 

1 1 X 

' 4 ’ 2 ’ 16 ’ 16 


2i H 3_]^/2 
2f H 3®3/2 
2f Kl 3f/2 
2® Kl 3®5/2 
2i B 83/2 
2i H 3_7/2 
2i Kl 85/2 
2i Kl 3^/2 


04/5 

fi® 

0-4/5 

fi® 

O16/5 


6 


-16/5 


6 


3/10 


6 


-3/10 
67710 
6 ® 


-■-7/10 

s® 

-’13/10 

B 

-13/10 

s® 

■’17/10 


fi® 

O-17/10 


6 : 


23/10 


-23/10 


6 : 


27/10 


-27/10 


6 : 


33/10 


-33/10 


6 : 


37/10 


-37/10 


1.7200 

1.7200 

1.7200 

1.7200 

1.9276 

1.9276 

1.9276 

1.9276 

1.9276 

1.9276 

1.9276 

1.9276 

1.9276 

1.9276 

1.9276 

1.9276 

1.9276 

1.9276 

1.9276 

1.9276 


10.854 

10.854 

10.854 

10.854 

14.472 

14.472 

14.472 

14.472 

14.472 

14.472 

14.472 

14.472 

14.472 

14.472 

14.472 

14.472 

14.472 

14.472 

14.472 

14.472 


1.1.1, C3\C3\C3 

1 . 1 . 1 , Cl, Cl, Co 
1 , 1 , 1 , Cl, Cl, Cl 

1.1.1, cl, cl, cl 

1.1, cl, cl, cl, del 
1,1, cl, cl, cl, del 
1,1, cl, cl, cl, del 

1.1, cl, cl, cl, del 

1.1, cl,d,cl,clcl 

1.1, cl, cl, cl, del 

1.1, ci,d,cl,clcl 

1.1, cl, cl, cl, del 
1,1, cl, cl, cl, del 
1,1, cl, cl, cl, del 

1.1, cl, cl, cl, del 

1.1, cl,d,cl,clcl 

1.1, cl, cl, cl, del 

1.1, cl,d,cl,clcl 

1.1, cl, cl, cl, del 

1.1, cl,d,cl,clcl 


_ 1 
3’ ■ 
1 1 
3 ’ 3 ’ 


1 J_ 2 

■ 3 ’ 15’ 15’ 5 

_J__1_ _2 

15’ 15’ 5 


2i4/5 ® 3-2 


^-14/5 


H 3, 


_ 1 _ 1 
3’ 3’15’ 

1 1 _ 4 _ _ 

3’ 3’ 15 ’ 

15 1 


A _2 

15 ’ 5 

J_ 2 
15 ’ 5 

2 7 


2^34/5 H 3£2 

2i4/5 ^ ^2 


2 ’ 16 ’ 10 ’ 5 ’ 80 

1 J>_ J_ _2_ 1_ 

2 ’ 16 ’ 10 ’ 5 ’ 80 

1 JL A _2 _3_ 

2 ’ 16 ’ 10 ’ 5 ’ 80 


^14/5 


KI3 


-5/2 


^-14/5 

1® 

^-14/5 


KI3, 
Kl 3. 


'5/2 

B 

7/2 


1 

7 


1 

2 


2 ’ 

16 ’ 

10 ’ 

5 ’ 


1 

3 


1 

2 

17 

2 ’ 

16 ’ 

10’ 

5 ’ 

80 

1 

3 1 


2 


17 

2 ’ 

16 ’ 10 ’ 

5 

, 

80 

1 

7 

1 


2 

13 

2 ’ 

16 ’ 

10 

, 

5 ’ 

80 

1 

7 

1 

2 


13 

2 ’ 

16 ’ 

10 

’ 5 

, 

80 

1 

1 


1 

2 

27 

2 ’ 

16 ’ 


10’ 

5 ’ 

80 

1 

1 1 


2 


27 

2 ’ 

16’ IQ 

1 ’ 

5 

, 

80 

1 

5 

1 


2 

23 

2 ’ 

16 ’ 

10 

, 

5 ’ 

80 

1 

5 

1 

2 


23 

2 ’ 

16 ’ 

10 

’ 5 

, 

80 

1 

1 

1 

2 

37 

2 ’ 

16 ’ 

10 

’ 5 

’ 80 

1 

1 

1 


2 


2 ’ 

16 ’ 

10 

, 

5 ’ 

f 

1 

3 

1 


2 

33 

2 ’ 

16 ’ 

10 

, 

5 ’ 

80 

1 

3 

1 

2 


33 

2 ’ 

16 ’ 

10 

’ 5 

, 

80 


2f4/5®3f,/3 


^14/5 


KI3 


-3/2 


M 9 

on " 


^-14/5 

B 

-14/5 


^3. 


'3/2 


^3 


-712 


2f4/5^3f/2 


^14/5 


^3 


- 1/2 


9 

on " 


^-14/5 

B 

-14/5 

1® 

=^14/5 


EI3 


1/2 


KI3 

KI3, 

KI3f 


-5/2 

B 

5/2 


^14/5 ^^1/2 

2-14/5 ^ 


•A 9 

on 


-14/5 

1® 

’^14/5 


KI3 

KI3 


-3/2 

■3^/2 


6 


fi® 

O1/7 

fi® 

0-1/7 

fi® 

O15/7 

B 

-15/7 


2.1082 

2.1082 

2.1082 

2.1082 


18.591 

18.591 

18.591 

18.591 


1.1, C5,C5,C5",C5 
I5 I5 CS; C55 C55 cl 

1 . 1 , Cl, Cl, Cl, Cl 

1 , 1 , cl, cl, cl, cl 


_i _i _n j_ 2 
4’ 7’ 28’ 28’ 7 

1 1 n__2 

4 ’ 7 ’ 28’ 28’ 7 

1 ^ _1 2 
4’ 28’ 7’ 7’ 28 

_1 __3_ 1 _2 13 

4’ 28’ 7’ 7’ 28 


2^1 
2f M 

3f/7 


2f S3f 


2^1 K1 3 


B 
-8/7 


'^B ,a 

30 

nB,b 

60 

o .B . Qr 

34 

RB,b 


2.1609 

2.1609 

2.1609 

2.1609 


20 

20 

20 

20 


1.1.2.2, V5,d5 

1.1.2.2, d5,d5 

1.1.2.2, d5,d5 

1.1.2.2, V5,d5 


1 _i 0 1 

5 ’ 5 ’ 2 

1 _1 1 _1 
5 ’ 5 ’ 4 ’ 4 

2 _2 Q 1 

5 ’ 5 ’ 2 


2 1 
' 5 ’ 4 ’ ' 


primitive 

primitive 

primitive 

primitive 


'^58/35 

^-58/35 

-B 

-’138/35 

B 

-138/35 


2.5359 

2.5359 

2.5359 

2.5359 


33.632 

33.632 

33.632 

33.632 


1, cl, cl, cl, cl cl, cl cl 
i,cl,cl,cl,clcl,clcl 
1, cl, cl, cl, cl cl, cl cl 
i,cl,cl,cl,clcl,clcl 


2 1 
5 ’ 7 ’ 
_ 2 
5’ ■ 

2 

5 ’ 


2 ^ 

7’ 35’ 35 

_i 2 m 
7’ 7’ 35’ 35 

1 2 _9_ _n 
7’ 7’ 35’ 35 

2 1 _2 n 9 

5’ 7’ 7’ 35’ 35 ^ 


^14/5 


^3 


-8/7 


— 2 
oc 


-14/5 


S3, 


'8/7 


2i4/5 B Sgyy 

-14/5 ^ gyy 


6 : 


46/13 


-46/13 


2.9132 

2.9132 


56.746 

56.746 


1, Cii, Cii, Cii, Cii, Cii 
1, Cii, Cii, Cii, Cii, Cii 


2 ^ 

’13’ 

4 

13 ’ 


-A A 

13 ’ 13’ 

_ 2 _ A _ 

13 ’ 13’ 


_ J_ 

13 

A A 

13 ’ 13 


primitive 

primitive 


6 


8/3 

3 

-8/3 


3.1107 

3.1107 


74.617 

74.617 


1 , C16, C16, C16, C16, C16 
1 , C16, C16, C16, C16, C16 


7 3 + V ^ S-TvAT O+yAT 

2 ’ 2 ’ 2 ’ 

1 3 + V ^ 3+V^ 3 + V ^ 

2 ’ 2 ’ 2 ’ 


5+V^ 7+V^ 
2 ’ 2 
5 + V ^ 7 + V ^ 


1 

’ 9 ’ 
1 

■ 9 ’ 


1 1 
9’ 3’ 
_1 
9 ’ 


primitive 

primitive 


6? 

6^2 


3.3263 

3.3263 


100.61 

100.61 


1 _2 3 n 1 

7’ 7 ’ 7’ ’ 3 

3 n 1 


7 

1 2 
7 ’ 


- -- 0 -- 
7, 7, 3 


primitive 

primitive 















































15 


TABLE III: A list of all 24 bosonic rank N = 7 topological orders in 2+ID with max(A^^'’) < 1. Since A = 7 is a prime 
number, all those 24 topological order are primitive. 


Af 

‘S'top 


di, d 2 , •■■ 

Si, S 2 , • • • 

’jB,a 
' 2 

1.4036 

7 

1 , 1 , 1 , 1 , 1 , 1,1 

07 

12 2 
75757 

5 

3 

7’ 

3 

7 

rjB, a 
^ -2 

1.4036 

7 

1 , 1 , 1 , 1 , 1 , 1,1 

0 ,- 

1 1 

75 75 

2 

7 


2 3 3 

r’ 7’ 7 

7B 
' 1/4 

2.3857 

27.313 

li IjCeiCeiCliCliCl 

0 ,| 

5 

32 ’ 

5 

32 

1 

4 ’ 

1 7 

4’ 32 

7B 

'-1/4 

2.3857 

27.313 

li liCeiCeiCeAeAl 

0,1 

5 5 

32 ’ 32 ’ 

1 

4 ’ 

1 

4 

7 

32 

7B 
' 3/4 

2.3857 

27.313 

li liCeiCeiCeAeAl 

0,1 

9 9 

32 ’ 32 ’ 

1 

4 ’ 

1 

4 

3 

32 

7B 

'-3/4 

2.3857 

27.313 

li liCeiCeiCliCliCl 

0,1 

9 

32 ’ 

9 

32 

1 

4 ’ 

1 3 

4 ’ 32 

7B 
' 5/4 

2.3857 

27.313 

li liCeiCeiCliCliCl 

0,1 

1 

32 ’ 

1 

32 

1 

4 ’ 

1 11 

4’ 32 

7B 

'-5/4 

2.3857 

27.313 

li liCeiCeiCliCliCl 

0,1 

1 1 

32 ’ 32’ 

1 

4 ’ 

1 

4 

11 

32 

7B 
' 7/4 

2.3857 

27.313 

li liCeiCeiCeAeAl 

0,1 

13 13 
32 ’ 32’ 

1 

4 ’ 

1 

4 

1 

32 

7B 

'-7/4 

2.3857 

27.313 

li liCeiCeiCeAeAl 

0,1 

13 

32 ’ 

13 

32 

1 

4 ’ 

1 1 

4 ’ 32 

7B 
' 9/4 

2.3857 

27.313 

li liCeiCeiCliCliCl 

0,1 

3 3 

32 ’ 32’ 

1 

4 ’ 

1 

4 

15 

32 

7B 

'-9/4 

2.3857 

27.313 

li liCeiCeiCliCliCl 

0,1 

3 

32 ’ 

3 

32 

1 

4 ’ 

1 15 

4 ’ 32 

7B 
' 11/4 

2.3857 

27.313 

li liCeiCeiCeAeAl 

0,1 

15 

32 ’ 

15 

32 

1 

4 ’ 

1 5 

4 ’ 32 

7B 

' -11/4 

2.3857 

27.313 

li liCeiCeiCeAeAl 

0,1 

15 15 
32 ’ 32’ 

1 

4 ’ 

1 

4 

5 

32 

7B 
' 13/4 

2.3857 

27.313 

li liCeiCeiCliCliCl 

0,1 

7 7 

32 ’ 32’ 

1 

4 ’ 

1 

4 

13 

32 

7B 

'-13/4 

2.3857 

27.313 

li liCeiCeiCliCliCl 

0,1 

7 

32 ’ 

7 

32 

1 

4 ’ 

1 13 

4’ 32 

7B 
' 15/4 

2.3857 

27.313 

li liCeiCeiCeAeAl 

0,1 

11 

32 ’ 

11 

32 

1 

4 ’ 

1 9 

4 ’ 32 

7B 

'-15/4 

2.3857 

27.313 

li liCeiCeiCeAeAl 

0,1 

11 11 
32 ’ 32’ 

1 

4 ’ 

1 

4 

9 

32 

rjB,h 
' 2 

2.4036 

28 

1,1,2,2,2,77,77 

0 , 0 , 

1 2 
7575 

3 

7’ 

L 

B ’ 

3 

8 

yB ,c 
' 2 

2.4036 

28 

1 , 1 , 2 , 2 , 2 , 77,77 

0 , 0 , 

1 2 
7575 

3 

7’ 

1 

8 ’ 

3 

8 

ijB ,c 
'-2 

2.4036 

28 

1 , 1 , 2 , 2 , 2 , 77,77 

0 , 0 , 

1 2 
75 

3 
’ 7 

1 

’ 8 ’ 

3 

8 

rjB^b 
' -2 

2.4036 

28 

1 , 1 , 2 , 2 , 2 , 77,77 

0,0, 

1 2 
7’ 7 

3 
’ 7 

5 

1 3 

8 ’ 8 

7B 

'8/5 

3.2194 

86.750 

1 1 C 13 ! C 13 :Cl3 I C 13 : C 13 1 C 13 

0,- 

1 A 0 

2 

5 ’ 

1 

3’ 

1 

5 

7B 

'-8/5 

3.2194 

86.750 

1 1 Cw ! C 131 C 13 I C 13 1 C 13 1 C 13 

o.| 

— — 0 
15’ ^ 


2 

5 ’ 

1 1 
■ 3 ’ 5 


topological order to an topological order ,Si,c) only 
shift the central charge c by 8: {Nl\si,c) —)■ {Nl^ ,Si,c + 
8). Such an operation is invertible. 

In our lists of 2+ID topological orders, we will only 
list topological orders up to invertible topological orders, 
i.e. we will only list the quotient 

{Topological ordersj/jlnvertible topological orders}. 

(48) 

It turns out that modular tensor category only describe 
topological orders up to invertible topological orders. 


C. A list of 2+lD bosonic topological orders with 
rank A = 1, 2, ■ • ■ ,7 

Table I lists 2+ID bosonic topological orders with rank 
N = 1,2, • • • ,5 and with = 0,1, 2, 3. Here we have 
ignored the invertible topological orders.®^ So the term 
“topological order” really refers to topological order up 
to invertible topological orders. 

In the table, there is 1 rank A = 1 topological order, 
which is actually a trivial topological order (i.e. corre¬ 
sponds to many-body states with no topological order). 


There are 4 non-trivial rank N = 2 topological orders, 
which correspond to = 1/2 bosonic Laughlin state with 
central charge c = 1 and the Fibonacci state with central 
charge c = ^, plus their time reversal conjugates. Those 
4 topological orders orders are primitive in the sense that 
they cannot be obtained by stacking non-invertible topo¬ 
logical orders with lower rank. 

Our numeric calculation also produce 12 rank N = 

3 and 10 rank A = 5 topological orders, which are all 
primitive since A = 3, 5 are prime numbers. 

For rank A = 4 topological orders, we find 18 of them. 
Applying eqn. (46), we hnd that by stacking two of the 
rank A = 2 topological orders, we can obtain 3 -I- 3 -I- 

4 = 10 distinct rank A = 4 topological orders. (If two 
(A/^,Si,c)’s are the same up to a permutation of the 
indices, we will say they describe the same topological 
order.) Indeed, 10 of 18 rank A = 4 topological orders 
are not primitive, corresponding to the stacking two of 
the rank A = 2 topological orders (see the blue entries in 
Table I). We also see that 6 primitive topological orders 
are Abelian since their topological excitations all have 
unit quantum dimensions di = 1. There are only two 
non-Abelian rank A = 4 topological orders, which are 
related by time reversal transformation. 

We like to pointed out the Ref. 103 gives a complete 
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classification of all 70 modular tensor categories with 
rank TV < 4. Compare with such a classification result, 
we find that our list for 35 rank < 4 topological orders 
is complete. (The other 35 modular tensor categories 
have iSii < 0 and do not correspond to unitary theory.) 

We find 50 rank N = 6 topological orders with < 2 
(see Table II). Most of those 50 topological orders are not 
primitive and can be obtained by stacking rank N = 2 
and rank N' = 3 topological orders (see the last column 
of Table II), where we have denoted the topological orders 
by their rank N and their central charge c: ). Only 10 

among the 50 are primitive. We also find 24 rank N = 7 
topological orders with = 0,1 (see Table III). They 
are all primitive since 7 is a prime number. 


TABLE IV: The fusion rule j ®i for topological order = 
64 ^, which is same as the fusion rule of 50(5)2 current algebra. 
For example, a(g)a = l©a©/3. 


Si 

di 

0 

1 

0 

1 

2 

5 

2 

2 

5 

2 

0 

V5 

1 

2 

V5 

j\l 

1 

a 

a 

p 

7 

X 

1 

1 

a 

a 

p 

7 

X 

a 

a 

1 

a 

p 

X 

7 

Q 

a 

a 

l©a©/3 

a®P 

7® X 

7® X 

p 

p 

p 

a©/3 

1 © a © a 

7® X 

7® X 

7 

7 

X 

7® X 

7® X 

l©a©/3 

a © a © /3 

X 

X 

7 

7® X 

7® X 

a © a © /3 

l©a©/3 


D. Understand the topological orders in the lists 

1. Non-Abelian type of topological order 

In this section, we like to gain a better understanding 
of the topological orders in the lists. Let us first use 
the stacking operation to introduce the notion of non- 
Abelian type of topological order. Two topological order 
Cl and C 2 have the same non-Abelian type iff there exist 
Abelian topological orders Ai and A 2 such that 

Cl K „4i = C 2 Kl A 2 - (49) 

The quantum dimensions in Abelian topological orders 
are all equal to I, so topological orders with the same 
non-Abelian type must have the same spectrum of the 
quantum dimensions (disregard the degeneracy). 


We like to remark that eqn. (29) can be rewritten as 
= (52) 

n 

Since Ni commute with each other, their largest positive 
eigenvalues di satisfy 

djdk = E (53) 

n 

Thus, if we express the quantum dimension di in the ba¬ 
sis of algebraic numbers, such as C™, with integer coeffi¬ 
cients, we can see the fusion rule from the product 
of dfs. 


3. Topological orders of parafermion non-Abelian type 


2. Quantum dimensions as algebraic numbers 

We next note that the quantum dimensions are alge¬ 
braic numbers (the roots of polynomial with integer co¬ 
efficients), since they are eigenvalues of integer matrices. 
So it is helpful to express those quantum dimensions in 
terms of algebraic expressions, such as y/n. But -y/n is 
not enough. So here we introduce another set of algebraic 
numbers 


sin[7r(TO + I)/(n + 2)] 
sin[7r/ {n + 2)] 


(50) 


It turns out that we can express all the quantum dimen¬ 
sions that we find in terms of C™ and ^/n. 

We note that the quantum dimensions that appear in 
Z„-parafermion CFT theory^^^ are all given by C™. Also, 
the Z„-parafermion theory has a central charge 




2n — 2 
n + 2 


(51) 


This suggests that many topological orders that we ob¬ 
tain are related to Z„-parafermion theories. 


Using the above concepts, we see that that the two 
N = 2 non-Ableian topological orders have the non- 
Abelian type of the Z 3 -parafermion theory since their 
quantum dimensions contain (^ 3 . Similarly, the TV = 3 
topological orders have non-Abelian types of the Z 2 and 
Zs-parafermion theories. The primitive TV = 4 non- 
Abelian topological order has a non-Abelian type of the 
Zy-parafermion theory. Among the N = 5 topological 
orders, we see the non-Abelian types of the Z 4 - and Zg- 
parafermion theories. Among the primitive N = 6 topo¬ 
logical orders, we see the non-Abelian types of the Zn- 
parafermion theories. For N = 7, we see that there are 
16 topological orders with the non-Ableian type of the 
Zg-parafermion theory. 


4- Topological orders of SO{k )2 non-Abelian type 

However, there are four TV = 6 topological orders (see 
Table IV) and four N = 7 topological orders that are 
not related to the parafermion theories. They are the so 
called TY{A,x,t)^^ category studied in Ref. 114, with 
A = Z 5 for TV = 6 cases and A = Zy for TV = 7 cases. 
They belong to metaplectic modular categories, which are 
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TABLE V: The fusion rule j for topological order 


Si 

di 

0 

1 

1 

7 

Cl 

1 

7 

Cl 

1 

7 

C?2 

3 

7 

C!2 

j\i 

1 

Q 

P 

7 

X 

1 

1 

a 

P 

7 

X 

a 

a 

/3©7 

1©X 

P®x 

a © 7 © X 

p 

p 

1©X 

a © 7 

a © X 

d © 7 © X 

7 

7 

P®X 

a © X 

1 © 7 © X 

a © d © 7 © X 

X 

X 

a © 7 © X 

/3 © 7 © X a 

© /3 © 7 © X 

l©a©/3©7©2x 

TABLE VI: The fusion rule j ® 

i for topological order A® = 6^/3 


Si 

di 

0 

1 

9 

Cl 

1 

9 

Cl 

9 

Cl 

1 

3 

CiV 

1 

3 

Cle 

jV 

1 

a 

P 

7 

X 

V 

1 

1 

a 

P 

7 

X 

V 

a 

a 

1 © a © X 

7 © 7 

/3 © 7 

a © X © 7 

/3 © 7 © X © 7 

p 

p 

7 © ^ 

1©^©X 

a © ry 

P®X®V 

a © 7 © X © 7 

7 

7 

p®v 

a(Sri 

1 © 7 © X 

7 © X © 7 

a © d © X © 7 

X 

X 

a © X © 7 

/3 © X © 7 

7 © X © 7 

l©a©d©7©X©7 

a©d©7©X©2J7 

V 

V 

/3 © 7 © X © 7 

a © 7 © X © 7 

a © /3 © X © >7 

a©d©7©X©2r7 

l©a©/3©7©2x©2r; 


defined as any modular category with the same fusion 
rules as SO{k )2 for k odd. They have rank N = (fc + 7)/2 
and dimension = 4A^. They have two 1-dimensional 
objects and two -yn-dimensional objects objects. The 
remaining objects have dimension 2.^^® We also like 
to point out that the four N = 6 topological orders and 
the four N = 7 topological orders that are closely related 
to U{l)k/Z 2 orbifold CFT with k = 5,7.“®’“^ 


A. Abelian topological orders 

All the Abelian topological orders can be describe by 
the A-matrix and can be realized by multilayer FQH 
states. 

1. The topological order = 2f in Table 1, is de¬ 
scribed by a 1-by-l AT-matrix K = (2). It re¬ 
alized by the Laughlin wave function for bosons 


5. Other topological orders beyond parafermion non-Abelian 

type 


In addition to the SO{k )2 non-Abelian topological or¬ 
ders, there are also a few topological orders that are be¬ 
yond parafermion non-Abelian type (see Tables V and 
VI). Some of the fusion coefficient Nl^ = 2 for those 
topological orders. 


VI. PHYSICAL REALIZATION OF THE 
TOPOLOGICALLY ORDERED STATES 


2. The topological order 4f is described by another 
1-by-l AT-matrix K = (4), and is realized by 
the Laughlin wave function for ~ 


3. The 3® topological order is described by a 2-by- 


2 AT-matrix AT = 


2 1 
1 2 


and can be realized by 


a double-layer bosonic FQH state T'gs = - 


4. Stacking two 3® topological orders give rise to a 
topological order described by 


In this section, we will discuss some physical realiza¬ 
tion of the topological orders that we find through the 
classifying theory. In this section, we will refer different 
topological orders by their rank N and central charge c, 
and use to denote them. 


/2 1 0 0 \ 
12 0 0 
0 0 2 1 
Vo 0 1 2/ 


( 54 ) 
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Such a topological order has 9 different types of 
topological excitations. Their spins are given by 

{si} = { 0 , - X 4 , —- X 4 , }. ( 55 ) 


i.e. there are 4 types of topological excitations with 
spin and 4 types of topological excitations with 
spin — 


5. There are two Abelian 4^ topological orders. The 
first one is the Z 2 topological order described by 

K = ^2 realized by Z 2 spin 

liquids^®’®^ or toric code model.The other is 
the double-semion topological order described by 

AT = , which can be realized by a string- 

net model^'. 


6. The 5q topological order is described by AT = 
, and can be realized by a double-layer 
bosonic FQH state = n(^i - - 



7. The Abelian 7f topological order in Table III is 


described by AT = 


4 3 
3 4 


and can be realized by 


a double-layer bosonic FQH state = Y\{zi — 

8. The Af and 5f topological orders are described by 


/2 1 1 1\ 


/2 1 1 1\ 

12 0 0 

K^B ~ 

12 10 

10 2 0 


112 0 

Vl 0 0 2yi 


i^l 0 0 2/ 


( 56 ) 


They can be realized by a four-layer FQH states. 


B. Non-Abelian topological orders of 
Z„-parafermion type 

Most non-Abelian topological orders that we found are 
of the Z„-parafermion^^^ type. For such kind of Z„- 
parafermion-type non-Abelian topological orders all the 
quantum dimensions are of the form for a set of m’s. 
(Note that the quantum dimensions can be Cn = !•) 
this section, we will discuss the physical realization of 
some of the Z„-parafermion-type non-Abelian topologi¬ 
cal orders. 

1. The 3^^2 topological order in Table I is of the Z 2 - 
parafermion type (see Table VH) . It can be realized 


TABLE VII: The fusion rule ji (g) i for a Z 2 parafermion topo¬ 
logical order — 3^2- Such a topological order can be 
realized®®’®'* by wave function = [i’ 2 ^{{zi})]^. The 

edge states are described by SU{ 2)2 x 17(1) Kac-Moody alge¬ 
bra. Note that = \/2. 


Si 

di 

0 1 A 
^2 16 

1 1 Ci 

j\i 

1 'i/i a 

1 

1 'i/i a 


ij} 1 a 

a 

a a 1 Q) tp 


TABLE VIII: The fusion rule j ® i for a Z 3 parafermion 
(Fibonacci) topological order ~ ^ 21 / 5 - Such 

a topological order can be realized®®’®* by wave function 
4^4^ = ■ The edge states are described by 

517 ( 3)2 X 17(1) Kac-Moody algebra. 


Si 

di 

0 

1 

1 

4 

1 

7 

20 

cl 

2 

5 

Cs 

j\i 

1 

a 

a 

T 

1 

1 

a 

a 

T 

a 

a 

1 

r 

a 

a 

<7 

r 

1 ©T 

a(B a 

T 

T 

a 

a © cr 

l©r 


by the following filling-fraction v = \ bosonic FQH 
wave function whose non-Abelian properties was 
first revealed in Ref. 30,34 (Feb. 1991): 

(57) 

where 'I'J^*"({zj}) is the fermionic wave function 
of n filled Landau levels. The 'I'os state was 

■^5/2 

shown to be a non-Abelian FQH state described by 
517 ( 2)2 X 17(1) Kac-Moody current algebra, which 
is the same as Z 2 -parafermion xl7(l) x 17(1) non- 
Abelian FQH state. Ref. 30,34,119 also studied the 
fermionic version of the above Z 2 -parafermion non- 
Abelian state 

(58) 

with rank N = 6, central charge c = 5/2, and 
filling-fraction 1 / = 1/2. 


TABLE IX: The fusion rule j ® i for the simplest Z 3 
parafermion (Fibonacci) topological order = 2 (^/ 5 . 



0 f 

di 

1 cl 

j\-i 

1 a 

1 

1 a 

a 

a 1 © a 
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TABLE X: The fusion rule j ® i for topological order 
~ It can be realized®'^’^'^ by wave func¬ 
tion The edge states are described by 

^55/7 

517 ( 5)2 X 17(1) Kac-Moody algebra. 


Si 

di 

0 

1 

1 

4 

1 

1 11 

7 28 

C 5 Cs 

1 

28 

Cl 

2 

7 

Cl 

iV 

1 

a 

/ 

a 7j 

r 

/ 

T 

1 

1 

a 

/ 

77 77 

r 

/ 

T 

a 

a 

1 

/ 

77 77 

/ 

T 

T 

a 

a 

/ 

77 

1 © r' a © T 

cr' © T 

cr © t' 

/ 

a 

/ 

a 

77 

a © T 1 © t' 

cr © t' 

cr' © T 

T 

T 

/ 

T 

ct' © r (7 © r' 

1 © cr © t' 

0 © cr' © T 

/ 

T 

/ 

r 

T 

(T © r' cr' © r 

a © cr' © r 

1 © cr © r' 


2. The 3^y2 topological order in Table I is also of the 
Z 2 -parafermion type. It can be realized by the 
following filling-fraction = 1 bosonic FQH wave 
function 


^-2 


■^3/2 


= Al(- 


Zl - 2:2 -23 - -2:3 


\Yliz,-Zj)e 


-lEh.r 


(59) 


It is closely related to the tz = 1/2 fermionic 
Pfaffient state 63/2 first proposed in Ref. 35 (Aug. 
1991), which has a rank N = 6 and a central charge 


3. 
2 ■ 




3/2 


= Al( 


zi - Z2Z3- Z3 


'n( 


Zi-Zjfe 


(60) 


The above two Z 2 parafermion states (one for 
bosonic electrons and one for fermionic electrons) 
can also be described by patterns of zeros (or ID 
occupation patterns)®^^®^ {n/} = {noi rii, 77 . 2 , • • • }: 

vI/ 3 .^ :{n 7 } = 20 | 20 | 20 |-- - , 

^ 63 /, -{ni} = 1100 | 1100 | 1100 | • • • . (61) 


3. The ~ ^ 21/5 topological order in Table I 

is of the Zs-parafermion (or Fibonacci) type (see 
Table VIII). It has the same Z 3 -parafermion (Fi¬ 
bonacci) non-Abelian type as the topologi¬ 

cal order (see Table IX). The topological or¬ 
der can be realized by the following filling-fraction 
12 = 3/2 bosonic FQH wave function with non- 
Abelian properties:®°’®^ 

= [4^^"{{a})]^. (62) 

The tk/B state was shown to be a non-Abelian 

21/5 

FQH state whose edge excitations are described by 
517 ( 3)2 X C7(l) Kac-Moody current algebra with 


central charge c = ^. 30 , 34,119 level-rank 

duality, the 517(3)2 non-Abelian type is the same 
as the 577 ( 2)3 non-Abelian type, which is also the 
same as the Z 3 -parafermion non-Abelian type. The 
fermionic version of the above Z 3 -parafermion non- 
Abelian state is given by^^’^'^bis 

Vk = vl//L({^J)[VI7LL({^J)]2^ (63) 

which has rank N = 10, central charge c = 21/5, 
and filling-fraction iz = 3/5. Ref. 34 also con¬ 
structed/studied those type of non-Abelian FQH 
states using parafermion CFTs in 1992. The non- 
Abelian excitations from such non-Abelian FQH 
states can perform universal topological quantum 
computations. 

4. The 4^^g topological order in Table I is of the Z 3 - 
parafermion (Fibonacci) type. It can be realized by 
the following filling-fraction 72 = 3/2 bosonic FQH 
wave function described by the following pattern of 
zeros: {n/} = {no, «i, ^^ 2 , • ’ ’ 

vk 4 B : {ni} = 30|30|30| • • • , (64) 

9/5 

i.e. Ueven = 3 and nodd = 0. It is closely related to 
the 72 = 3/5 fermionic FQH state constructed using 
Z 3 parafermion CFT in 1998^^° with Nc = IO 9/5 
described by the pattern of zeros: 

^ 109/3 : = 11100|11100|11100| • • • . (65) 

5. The ^ topological order in Table H 

is of the Zs-parafermion type (see Table X). It 
can be realized by the following filling-fraction 
72 = 5/2 bosonic FQH wave function which is non- 

Abelian30’34. 

%,,, = ( 66 ) 

The fermionic version of the above Zo-parafermion 
non-Abelian state is given by®°’®^4i9 

vk = VI,LL({^J)[3l,LL({^^})]2^ (67) 

which has rank TV = 21, central charge c = 55/7, 
and filling-fraction 72 = 5/7. The above non- 
Abelian FQH states and their edge excitations are 
also described by 577 ( 5)2 x 77(1) Kac-Moody cur¬ 
rent algebra. 

6 . The 6 fgy.^ topological order in Table H is of the 
^ 5 -parafermion type. It can be realized by the fol¬ 
lowing filling-fraction 72 = 5/2 bosonic FQH wave 
function {ni} = {uq, ni, 772 , • • • }: 

:{’^/} = 50|50|50|--- . ( 68 ) 

It is closely related to the 72 = 5/7 fermionic Z 5 - 
parafermion state^^° with Nc = 2 I 15 / 7 : 

4'21i3/7 : = 1111100|1111100|1111100| • • • . (69) 
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C. Non-Abelian topological orders of 
Zn X 2„/-parafermion type 

Some non-Abelian topological orders that we found 
are of the Zn x Z„/-parafermion type. For such kind of 
Zn X Z„/-parafermion-type non-Abelian topological or¬ 
ders, all the quantum dimensions are of the form 
for a set of m, m'’s. Some of those topological orders can 
be realized by stacking Z„-parafermion topological order 
with Z„/-parafermion topological order. 

For example, staking two Za-parafermion topo¬ 

logical order described by wave function ^ Si^e 

us a third Z 3 x Za-parafermion ^ 2/5 topological 

order described by wave function 

{{zi},{Wi}) ='i/ 2 B ({z^})'^' 2 B ({wj). (70) 

28/5 14/5 14/5 


The spin spectrum has the {s^} —)■ {—Si} symmetry for 
all those topological orders. However, since c 7 ^ 0 im¬ 
plies a chiral edge state, those topological orders cannot 
be realized by time reversal symmetric systems. It was 
suggested in Ref. 121,122, that the topological 

order can be realized as the time-reversal symmetric sur¬ 
face states of a 3-l-lD time reversal symmetric symmetry- 
protected topological state. We believe all those topolog¬ 
ical orders can be realized as the time-reversal symmetric 
surface states of the same 3+lD time reversal symmetric 
symmetry-protected topological state. In other words, 
those topological orders have anomalous time-reversal 
symmetries, which have the same type of anomaly. 


F. 2-|-lD fermionic topological orders 


Similarly, staking Za-parafermion 2^^^^ topological or¬ 
der and Z 2 -parafermion 3^2 topological order together 
produce a third Z 3 x Z 2 -parafermion 6 ^ 3 ^ topological 
order in the Table II, which is described by wave function 


({Zi},{Wi}) = 4'2B ({2:j})4'3S ({rcj}). 

^33/10 14/5 ^1/2 ^ 


(71) 


We may identify Zi and Wi in the above wave function, 
trying to obtain a new topologically ordered state. If we 
are lucky, the new wave function 


(72) 

will describe a gapped state, which will be a topolog¬ 
ical order with one less central charge (for details, see 
Ref. 40), i.e. a Z 3 x Z 2 -parafermion 6 ^ 3^33 topological 
order. The Q 23 /W topological order does appear in our 
table II, which implies that identifying Zi and Wi will give 
us the Z 3 X Z 2 -parafermion topological order Q 23 /W 


D. 2-|-lD time-reversal symmetric topological 
orders 

We have found 6 topological orders with c = 0 and 
N < 7: three = 4^, one = 5^, and two 
= 6q . The spin spectrum has the {si} —>■ {—s^} 
symmetry for all those topological orders. It suggests 
that those topological orders can be realized by time re¬ 
versal symmetric systems. In contrast, the spin spectrum 
does not have the {s^} —)■ {—Si} symmetry for most topo¬ 
logical orders, suggesting that they cannot be realized by 
time reversal symmetric systems. 


E. 2-|-lD anomalous time-reversal symmetric 
topological orders 

We have found 4 topological orders with c = 4 and 
N <7: one = 4f, one = 5f, and two = 6 f. 


Although we have only discussed bosonic topological 
orders in this paper, we can see fermionic topological 
orders ^®424 q^j- classification of bosonic topological 
orders. Let us illustrate this point through an example. 

We start with the = 4®, Si = (0,0,0, \) topological 
order {i.e. the Z 2 topological order®®^®®) in Table I. We 
know that the Z 2 topological order contain a fermionic 
excitation /. If we add the fermionic excitations to the 
ground state and let the fermions to form a product state, 
such an addition will not change the Z 2 topological order. 
However, if we let the fermions to form ap-l- ip supercon¬ 
ducting state, then the Z 2 topological order will change 
to a different topological order. Since the p + \p super¬ 
conducting state has c = 1I2 edge state, the new topo¬ 
logical order should also has c = 1/2. This suggests that 
fermion condensation into the p + '\p state will change 
the = Aq Z 2 topological order to the iV^ = 3 (y 2 
topological order in Table I. We note that the Z 2 -charge 
and the Z 2 -vortex both behave like the same tt-Aux to 
the fermion /. In the p + ip state, tt-Hux will carry an 
Majorana zero mode and behave like a topological exci¬ 
tations of quantum dimension \/2 = / 2 - Such a kind of 
topological excitations appear in the = 3 f ^2 state, 
confirming our identification. 

Similarly, if we let the fermions to form 2n -|- 1 layers 
of p -I- ip superconducting states, then the Z 2 topological 
order will change to the A® = 3 ^_,_ 3^2 topological order 
in Table I. This is because 2n + l layers of p + ip states 
have chiral central charge c = ( 2 n -|- l )/2 edge state. 
Also, if we let the fermions to form 2n layers of p -I- ip 
superconducting states {i.e. a u = n integer quantum 
Hall state), then the Z 2 topological order will change to 
the = 4^ topological order in Table I. This is because 
2n layers of p -I- ip states have chiral central charge c = n 
edge state. 

We also see that the = 6 f( 3 +io„)/io states are 
related by the fermion condensation into v = An integer 
quantum Hall state. The = 7 ^( 3 _,_ 2„)/4 states are 
related by the fermion condensation into An layers of 
p -I- ip states. 
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VII. A CLASSIFICATION OF 1+lD 
GRAVITATIONAL ANOMALIES 

Since the 1+lD bosonic gravitational anomalies (both 
perturbative and global gravitational anomalies of known 
or unknown types) are classified by the 2+lD bosonic 
topological orders, {S, T, c) or ,Si, c) give us a classi¬ 
fication of all 1-l-lD bosonic gravitational anomalies. We 
may also view the tables I, II, and III as tables of simple 
bosonic gravitational anomalies. When c 7 ^ 0, the 1-l-lD 
bosonic gravitational anomaly contain perturbative grav¬ 
itational anomaly. When c = 0, the 1-l-lD bosonic gravi¬ 
tational anomaly is a pure global gravitational anomaly. 

Given a I-I-ID low energy effective theory £i_|_id, how 
do we know if the theory has gravitational anomaly or 
not? According to Ref. 63,123, we first try to realize 
£i+id by the edge of 2-l-ID gapped liquid system de¬ 
scribed by £ 2 + 10 - We then use the non-Abelian geo¬ 
metric phase^^ or wave function overlap®^ to compute 
5, T. From (S', T, c), we learn the type of the gravita¬ 
tional anomaly in the I-t-ID theory £i_|_id. 

As an example, let us consider the following I-I-ID 
bosonic system 

£ i - i-id = —;^dx(t>idt4>j — -^Vijdx4>idx4'J, ( 73 ) 
47r 47r 

where (pi are compact real fields {(pi ~ ())/ + 27r). Such 

1- l-lD effective theory can be realized by the edge of 

2- l-ID iF-matrix FQH state.We find that the I+ID 
effective theory £i+id is anomaly free if det(Ar) = ±I 
and K has an equal number of positive and negative 
eigenvalues. 

If K has different numbers of positive and negative 
eigenvalues, then the above I-I-ID bosonic theory will 

have a perturbative gravitational anomaly, li K = 

[2 0\ /2 3 \ 

, A = etc , then the above I-I-ID 

- 2 y \^3 2 y 

bosonic theory will only have a global gravitational 
anomaly. 


VIII. SUMMARY 

In this paper, we review the discovery and development 
of topological order - a new kind of order beyond Landau 
symmetry breaking theory in many-body systems. We 
stress that topological order can be defined/probed by 
measurable quantities {S,T,c) or {N{p, Si,c). 

We know that symmetry breaking orders can be de¬ 
scribed and classified by group theory. Using group the¬ 
ory, we can obtain a list of symmetry breaking orders, 
such as the 230 crystal orders in three dimensions. Sim¬ 
ilarly, in this paper, we present a simple theory of 2+ID 
bosonic topological order based on {S, T, c) or {Nl ^, Sz, c). 
This allows us to obtain a list of simple 2-l-ID bosonic 
topological orders. Although it is not clear if the theory 


presented in this paper is a complete theory for topolog¬ 
ical order or not, it serves as the first step in developing 
such a theory. 

We also discussed how to realize the some of the topo¬ 
logical orders in the list by concrete many-body wave 
functions. A more systematic way to realize those topo¬ 
logical orders is via simple current CFT, which will ap¬ 
pear elsewhere. 
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Appendix A: A fusion category theory for the 
amplitudes of planar string configurations 

In Section V, we simply list many conditions on 
{S, T, c) or {N]p, Si, c). We did not explain where do they 
come from, although they are derived in various math¬ 
ematical literature (for a review, see Ref. 104). In the 
next a few sections, we will try to explain and understand 
some of those conditions, in a simple and self-contained 
way. 

1. The string operators 

Although a topological excitation cannot be created 
alone, a pair of particle and anti-particle i,i can be cre¬ 
ated by an open string operator Wi. In some cases, the 
open string operator Wi is a product of local operators 
along the string 

vFz= n (^1) 

zGstring 

But more generally, the open string operator Wi has a 
more complicated structure. We need to use local op¬ 
erators with two “bond” indices, Mf^{xi) to construct 
it:'" 

W^= AI““^(U)M“'“^(z2)M“^“^(i3)--- (A2) 

o.ia.2a3--- 

We see that the bond indices are traced over and the 
above string operator is a matrix-product operator. 

We may choose the local operator M°'^{xi) properly 
such that the normal of Wi| ground) does not depend 
on the length of the string operator, and furthermore 
Wi I ground) cx [ground). Such a string operator obeys the 
so called “zero law” as described in Ref. 125. Ref. 125 
pointed out that it is always possible to obtain such “zero 
law” string operator for each type of topological excita¬ 
tion. 

Using the “zero law” closed string operator, we can 
have another way to understand the simple type and 
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FIG. 6: World-lines in a local region represent a local tun¬ 
neling process, where the topological excitations are created 
in pairs, and then braided and fused, and at last annihilated 
in pairs. The picture is also a 2D projection of a 3D string 
configuration. 

composite type. If z is of a composite type, i = j0fc0- • •, 
then the corresponding “zero law” closed string operator 
can be decomposed into a sum of “zero law” closed string 
operators: 

W, = W,+Wk + --- . (A3) 

If a “zero law” closed string operator cannot be decom¬ 
posed, then the corresponding particle is of a simple type. 
The correspondence between the composite type and the 
sum of the string operators, as well as the correspon¬ 
dence between the fusion of topological excitations and 
the product of string operators, allow us to see that the 
“zero law” closed string operators for simple types satisfy 
an algebra described by the fusion coefficients 

(A4) 

k 

(We will derive this relation later in Section A 12.) 

2. The space-time world lines and quasiparticle 
tunneling process 

In the space-time path integral picture, a “zero law” 
string operator correspond to a string in a time slice of 
a fixed time. We can consider more general “zero law” 
strings in space-time that can go through different times. 
Those more general strings in space-time correspond to 
the world-lines of the topological excitations. If the all 
the world-lines are confined in a local region, then they 
will represent a local tunneling process, where the topo¬ 
logical excitations are created in pairs, and then braided 
and fused, and at last annihilated in pairs (see Fig. 6). 
Since the degenerate ground state are locally indistin¬ 
guishable, such a local tunneling process causes the same 
amplitude for different degenerate ground states. Thus 
world-line confined in a local region correspond the a 
complex number (the amplitude) in the space-time path 
integral picture. Since the world-line correspond to “zero 


FIG. 7: A local tunneling process of two linked world-lines. 
The dash lines are the framing of the world-lines. In such a 
tunneling process the framing dash lines do not link with the 
original strings, indicating that the particles do not twist (or 
rotate) during the tunneling process. The amplitude of the 
above linked loops is a complex number denoted as 

law” strings, the above complex number (the amplitude) 
does not depend on the shape and length of the world¬ 
line. It only depend on the linking and the fusion of the 
world lines. 

Clearly, number of types of the strings is given by the 
rank N, and the strings are labeled by the simple type i of 
the topological excitations. Also, the strings are oriented 
if the particle i and anti-particle i are different (see Fig. 

The fusion of particles is represented by the branching 
point of the strings. If N'j^ = 0, then the amplitude of 
the string configuration that contain a branching point of 
i,j, k strings will be zero. When > I, it means that 
the space i® j contain several copies of the space k. We 
will label the copies of the space fc by a = 1,2, • • • , Nl^. 
There will be a tunneling amplitude into each copy of 
the space k, and the tunneling amplitude will depend 
on a. So we will include the index a € [1, • • • , Nl^] on 
each branching point. In this case, each such labeled 
graph of strings gives rise to an amplitude. We will use 
A{X) to represent such an amplitude for a labeled string 
configuration X, such as the one in Fig. 6. 

We also like to mention that the world-lines have a fi¬ 
nite cross section which is not circular. So more precisely, 
the world-lines are represented by framed strings in 7. 
The framing represents the finite cross section, which do 
not have rotation symmetry. 

3. Planar string configurations 

In this paper, we will mostly draw the string in 3- 
dimensional space-time in terms of their projection on a 
particular plane. In the 2D projected representation, we 
will always choose a canonical framing, by displacing all 
the 2D strings a little bit in a direction perpendicular to 
the 2D plane to obtain the framing dash-lines. So when 
we draw such 2D string configurations, we will assume the 
above canonical framing and will not draw the framing 
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dash-lines. 

In the rest of this section, we will consider all the ampli¬ 
tudes for planar string configurations, where the strings 
in the 2D projection do not cross each other. It turns out 
that the amplitudes for different planar string configura¬ 
tions have a lot of relations, so that we can determine the 
amplitudes for all the planar string configurations from a 
set a tensors that satisfy a certain relations. This turns 
out to be a fusion category theory of the amplitudes for 
planar string configurations. Ref. 78 presented such a 
fusion category theory for more general fermionic case. 
Here, we will present the simplified case for bosons. 

Since the planar string configurations can be viewed 
as the world line of particles in 1-l-lD space-time, the 
fusion category theory described here can also be viewed 
as the classifying theory for anomalous 1-l-lD topological 
orders,which can be described by the particles 
tunneling process in 1-l-lD space-time. 


4. The first type of linear relations: the F-move 


Let us consider a local region in the 2D projected 
string configuration. We fix all strings cutting across 
the boundary of the region, and consider all the differ¬ 
ent ways that the strings connect to each other in the 
region. Those different string configurations describe dif¬ 
ferent local tunneling process. If a subset of string config¬ 
urations already describe all the channel of the tunneling 
processes, then the amplitude of every other local string 
configuration can be expressed as a linear combination of 
the amplitudes for the subset of string configurations. 

‘ ^ 


In fact, the graph 


with fixed ijkl but different 


mafi is a subset of string configurations that describe all 
the channel of the local tunneling processes with fixed 

i j y^k 

ijkl. The graph \5^ with fixed fc/but different nyd 


is another subset of string configurations that describe all 
the channel of the local tunneling processes with fixed 
ijkl. So we can express the amplitudes for string con¬ 
figurations in one subset in terms of the amplitudes for 
string configurations in another subset: 


A j = i: FZ:Sa j . (A5) 

We note that 

Kin.xS^ = 0 when (A6) 

iV- < 1 or NU < 1 or < 1 or < 1. 

When < 1 or < 1, the left-hand-side of 

eqn. (A5) is always zero. Thus = 0 when < 1 


or < 1. When N-F < 1 or A”™ < 1, amplitude on 
the right-hand-side of eqn. (A5) is always zero. So we 
can choose = 0 when < 1 or A™ < 1. 

For fixed z, j, k, and I, the matrix F^j with matrix 
elements ^ matrix of dimension 

X The matrix describe the 

relation of the tunneling amplitude through one set of 
channels described by basis majj and through another 
set of channels described by basis nyd. We note that 
the tunneling maps i,j,k to I with degeneracy. The first 
tunneling path gives rise to basis maj3 of the degenerate 
subspace. The second tunneling path gives rise to basis 
nxS of the degenerate subspace. The degenerate sub¬ 
space of I should to the same, regardless the tunneling 
paths. So we require to satisfy 

= Y,NrNF^ ( a ?) 


and the matrix F^j to be unitary: 


^kin.^5 ^ (T'Jn.v/) — ^m,m' Sa,a'S . (A8) 

nxS 


(But here we do not require = A^*). It is easy to see 
that the unitary condition implies: 



Similarly, we have a dual F-move 



where also satisfies a unitary condition. 

The F-move (A5) can be viewed as a relationship be¬ 
tween amplitudes for different graphs that are only differ 
by a local transformation. Since we can transform one 
graph to another graph through different paths (i.e. dif¬ 
ferent sets of local F-moves), the F-move (A5) must sat¬ 
isfy certain self consistent conditions. For example the 

j k l i j 


graph 



can be transformed to 



through 


two different paths; one contains two steps of F-moves 
and another contains three steps of F-moves as described 
by eqn. (A5). The two paths lead to the following rela¬ 
tions between the wave functions: 
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(All) 


/i J k A 



qSe 



The consistence of the above two relations leads a condi¬ 
tion on the F-tensor: 

AT’* N*J 

pitn.c/SX 

/ ^ ^ ^ ^ kntj'qif ^ IpSjKj lsq,64> 

t r)—l if—l K=1 

N^c, 

= E F!;Sf’‘K"H- (A13) 

e=l 

which is the famous pentagon identity. The above pen¬ 
tagon identity (A13) is a set of nonlinear equations sat¬ 
isfied by the rank-10 tensor . The above con¬ 

sistency relations (A13) are equivalent to the require¬ 
ment that the local unitary transformations described 
by eqn. (A5) on different paths all commute with each 
other. 


and 

Qifc,a/3 ^ Q jy3k ^ ^ 

We note that the number of choices for the four indices 
(j, k, a, j3) in must be equal or greater than 1: 

A = )" ^ 1- (A17) 

jk 


6 . The third type of linear relations: the Y-move 


5. The second type of linear relations: the O-move 

The second type of linear relations re-express the am- 
1 / 

plitude for it terms of the amplitude for 


A 


^ X ^ 

/"a\ 

V ^ 


= of 


(AM) 


We call such local change of graph an O-move. Here 
03 k,oiP satisfies 

jsjjk jsjjk 

EEE 

kj a — 1 / 3—1 


/Wi 

For fixed i,j, f for different kap describe all the 
iXpXi 


possible tunneling channels. So the amplitude for 

can be expressed in terms of the amplitudes of f : 

VpXi 






= A 


k,oc^ 


(A18) 


We will call such a local change as a Y-move. We can 
choose 




(A19) 
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7. A relation between and 

We find that the following tunneling amplitude has two 
ways of reduction: 



(A20) 


M 


Pi 


Pi 



(A21) 


The two reductions should agree, which leads to the con¬ 
dition 


Qik,aX 


E vjk ^jk,j\^jk,aP 
^i,Pl'^i 

Pi 


(A22) 


8 . A freedom of changing basis at each vertex 

We note that the following transformation changes the 
basis at the branching point labeled by a 


The above transformation induce the following trans¬ 
formation on 


o 


,jk,al3 A,a pjk,l3 ^jk,a'/S' 

i ^ J jk,a.' •/ i,l3' ’ 

\ri3 , ( ( fk,^' 

^ k,a/3 w fc,a / ' '^k,a'l3' 


(A25) 


T^ijm.ot^ pij,ct 

^ kln,x^ ^ ^ rn^ot'J l,l3' \Jn,x ' \Jl,6 / ^kln,x"5' ' 


We note that the first line of the above equation is a 
singular value decomposition, since and for 

fixed i,j-,k, are independent unitary matrices. Tkus, we 
can use the above basis-changing freedom to choose 


qM ^ Q3k,<.. 


a.^1 


o 


\jk,a 


> 0 . 


(A26) 


We see that as the singular values, can be chosen 

to be positive real numbers. Then eqn. (A22) implies 
that 


= yor- 


(A27) 


9. A relation between and 


We also hnd another graph that can have two ways of 
reduction as well: 



_ pjkl,^iT A 

/ j mis,x<^ 


V / 



V / 


_ T^ikl,iiT Qkm,XQjP,0‘ A 
^ niip.xCA^P ^ 


(A28) 



= F: 


\jkl.,liT 


mip,x^ 


V / 



V / 


_ pjkl,pT fjjk,fj.Qlm,T A 

^ mip,xoi^l 


(A29) 


This allows us to obtain another condition 

pjkl,fiT _ pjkPpT Qkm,XQjP,af0lm,T-,-liQjk,fj,\-l 

^ mip,xoi. ^ mip,xoi^P ^ ' 

(A30) 


where fl^ is a unitary matrix 


We require to be unitary, which leads to 


Y.&(&'r=S.v. (A24) 

s 


Similarly, we have unitary transformation for ver¬ 
tices with two incoming edges and one outgoing edge. 
Such transformations correspond to a choice of basis and 
should be regarded as an equivalent relation. 


(jkm,x'Qjp',a' Qkm,x(l3P,°‘ 

E /pjkl,ii.T \*_fV_ 2 Tp3kl,p,T ^p k'i 

V mip',x'^ 03k,p mip^x^ 0lm,T^jk,^ 


E 

Ipr 


f pjkl,pT \ " 
\^mip' ,x'a') 


pjkl,pT 

mip,xoL 






(A31) 
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or 


E 


(F; 


<jkl,fiT \>ic 




Spp' ^aa' 




2 ’ 


(A32) 


The above condition can be satisfied by the following 
ansatz (note that is real and positive) 

where = 1 for > 0 and Sf^ = 0 for = 0. 
From eqn. (A15), we find that Wi satisfy 


E = WkD^, D = wf. (A34) 

v V ' 

The solution of such an equation gives us Wi. 

Let us consider the fusion of n type-i particles. The 
dimension of the fusion space is 


We see that cf" is the left eigenvector of Ni with eigen¬ 
value di'. 

d^N, = d^di. (A40) 

In other words, the left eigenvector of Ni with the largest- 
eigenvalue is independent of i. Such a common left eigen¬ 
vector is given by c?^, and the corresponding largest- 
eigenvalue is the quantum dimension di. 


10. The fourth type of linear relations: the H-move 


Let us consider a new type of move - iL-move. First, 


/'■ 


for fixed z, j, k, I 


lA. \ • 

' m \ 


A 


and 


V 


describe all 


the possible tunneling channels and their can express 
each other via unitary linear relations: 


D^{n)= Y (A35) 

Let di be the eigenvalue of matrix Ni (defined as {Ni)kj = 
Nl^) with largest absolute value, di will called quantum 
dimension of type-z particle. Since all the entry of Nl^ are 
non-negative, one can show that di is real and positive. 
We see that the dimension of the fusion space is given by 




a 




A 




V '' Tjkim,OL^ A 

is 

k/Y 
\ / 





nx5 


In the following, we will show how to compute the coef¬ 
ficients from and w,. 

First, by applying the Y-move, we have: 


i4*(n) =0(l)d”. (A36) 


Now, consider type-z particles and type-j particles. 
We first fuse n type-z particles, then fuse the result with 
n type-j particles, and then fuse the result with n type-z 
particles, etc . The dimension of the fusion space is 


E 

kl t" ' — 1 


Pj-ik^n — 2 J\^‘i'k2n J\^‘^k2n — 1 

k^n-l k2n+l k2n 


KT-T ■ ■ ■ (A37) 


The above dimension of the fusion space should be 
0{l)df d" . But if the largest-eigenvalue eigenvectors 

of Ni and Nj are different, we will get 0(l)ci”^d"^/” 
as the dimension of the fusion space. The fact that 
/ = 1 implies that the largest-eigenvalue eigenvectors 
of Ni and Aj must be the same (as implied by the condi¬ 
tion eqn. (14)). Let (z;i,U 2 , • • •) be the common largest- 
eigenvalue eigenvector for Aj’s: 


NiV = diV. 


(A38) 


Since all the entry of Aj are non-negative, one can show 
that Vi is real and positive. Using eqn. (14), we find 


^ NliNr'^Vk = Y 

m,k m,fe 


^ Nj^dmVi = ddjvu ddj = Y 


(A39) 


n 

ty 

A 


LA- 


\ 




E 


n,x'S 



(A42) 


Next, by applying an inverse F-move, we obtain: 





\ 


/ 

(A43) 


Finally, by applying the 0-move, we end up with: 



All together, we find: 

[ Tpkmi,P'x\^ I^km,ap' 

jnl,l3x' j 


H 


kim,OL(3 


jln,xS 


X'P' 


(A44) 


(A45) 
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Under the proper basis choice eqn. (A26), we can fur¬ 
ther express the coefficients as: 

T^kim,a(3 _ -Y'kl f T^kmi,ax\*^km,a 

^jln,xS ^n,S\^jnl,^S ) 

= (A46) 

The unitarity condition for H-move requires that: 


E 


T^km'i,CKX f Tpkmi,cxx 
jnl,0'5' V jnl,/36 




* 


1.' ^/1 




(A47) 


With the special ansatz eqn. (A33), we can further sim¬ 
plify the above two expressions as: 


h: 


kim,OL^ 


jln,xS 


WmWn 

WiWl 


( pkmi,ax\’’ 
\^jul,l3S ) 


(A48) 


and 


E^ 

nxS 


Tjkm'i,ax f T^kmi,ax\* r s- c 

i^jnl,/3S ) “ ~ Omm'Oaa'Opp 


(A49) 


Similarly, we can also construct the dual-H move: 


t 








k 


y'y 





(a) (b) (c) 


FIG. 8: (a): Two tunneling processes: Wf and W/. (b): 

The tunneling path of the above two tunneling processes can 
be deformed according to the Y-move (c): The O-move can 
reduce (b) to (c). 


Those conditions form a set of non-linear equations 
whose variables are ^ (where Wi can be 

determined by alone). Let us collect those conditions 
and list them below 




N 


N^^NV 


m—0 


N 

E 

n—0 



E: 



A 


m 
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and we can express x6 


as: 


frkim^Oi^ _ ^kim^a/S 

^jln,xS ~ ^jln,xS ’ 


(A50) 


(A51) 


where the coefficients 


can be expressed as: 


fTkim,a0 

^jln.xS 


_ Xrkl jpimk.aS^ml,l3 

— n,6^lnj,0x 


= F, 


imk,ocx (r^kl,5 


lnj,^6 






ml,^ 

j 


(A52) 


Again, with the special ansatz eqn. (A33), we have: 


. ^ 1; 

jk 

(A54) 


nxS 

• Kin,xS^ = 0 when 

Nii < 1 or <loT <1 or M;™ < 1, 
tv)'" Af“ Af“ 

^ pij‘m,a/3pitn,ipxpjkt,riK 

/ j / j / j / j knt,rnp IpSjKj ^ lsq,S<f) 
t rj—1 (p—1 K—1 

= E <A55) 

€=1 


^kim,aj3 _ ^imk,a6 

jiu,xs - Y WjWfe '"fAx 


(A53) 


It is easy to see that the unitarity condition for dual H- 
move is automatically satisfied if the H-move is unitary. 


• '^w^WjNl^ =WkD‘^, D= Y. ■ 

ij V * 


(A56) 


11. Summary of the conditions on the linear 
relations 




'km'i,a.x f ipli^^,0'Xy _ 


,OLX( -pk'i 

'5' V'^7U 


^ jnl,^'S' v'^ jnl,05 


nxS 


WiWl 


Wr, 


(A57) 


We see that valid tunneling amplitudes A(X) can 
be characterized by tensor data How¬ 
ever, only certain tensor data {Nk , that sat¬ 

isfy the conditions eqns. (A8, A6, A13, A34), can self- 
consistently describe valid tunneling amplitudes A{X). 


12. A derivation of string fusion algebra 

As an application of the above algebraic structure, let 
us consider the “zero-law” string operators on a torus 
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Si X Slj, wrapping around the 2 :-direction Wf. View¬ 
ing those string operators as world-lines in space-time, 
applying the Y-move and then the 0-move, and using 
eqn. (A27) (see Fig. 8), we find that 

w:w^ = 

k,a. 

= (A58) 

k 

We see that the algebra of the loop operator VFf forms 
a representation of fusion algebra i® j = N^^k. The 

operators Wf = SWfS~^, where S is the 90° rotation, 
satisfy the same fusion algebra 

WfWf = ^ Nl^f. (A59) 

k 

This way, we derived eqn. (A4). 


Appendix B: Unitary m-fusion category 

The tensor data satisfying the condi¬ 

tions eqns. (A54), (A55), (A56), and (A57) form a so call 
unitary m-fusion category. In fact, we can view the graph 

and the graph as the different ways to fu¬ 
sion three particle types i,j,k to = 

copies of particles 1. But the two ways of 
fusions lead to different basis of the space of copies 
of 1. The tensor is nothing but the unitary trans¬ 

formation that relates the two basis. 

However, here we do not require the existence of trivial 
particles type. Thus the structure we described is not a 
unitary fusion category ans we call it an unitary m-fusion 
category (UmFC). 

We like to stress that the fusion discussed here is not 
symmetric {i.e. we do not require N^f = Njf). Thus 
fusion that we are talking about is the fusion of ID par¬ 
ticles, where their order cannot be changed. Therefore 
UmFC is a classifying theory of 1-l-lD anomalous topo¬ 
logical orders Such anomalous topological or¬ 

ders cannot be realized by any well defined ID lattice 
models, but they can realized as boundary of 2D lattice 
models with non-trivial 2-l-lD topological orders. Those 
2-l-lD topological orders C 2+1 are described by modular 
tensor categories, which are uniquely determined by the 
1 -|- ID anomalous topological order on the boundary. 
In fact, the 2-l-lD bulk topological order is the Drinfeld 
center of the 1-l-lD anomalous boundary topological or¬ 
der: C 2+1 = Z{Ci+i). One concrete way to compute the 
Drinfeld center is described in Ref. 126. 


Appendix C: Unitary fusion category and the trivial 
particle type 

1. Trivial particle type and rule of adding trivial 
strings 

In the above discussion, we did not assume the exis¬ 
tence of a trivial particle type. Here we will assume the 
existence of such a trivial particle type, denoted by 1, 
which satisfies the following fusion rule 

l®i = i®l = i. (Cl) 

Thus N]f satisfies 

= (C2) 

We also requires that for every i there exists a unique i 
such that 

l = z, 1 = 1, (C3) 

By setting ^ = 1 in eqn. (14), we find the following sym¬ 

metry condition on N^f : 

NF = Nl’". (C4) 

Using the above, we can rewrite the condition eqn. (A56) 
as 

^ wiWjNl'^ =wj,D^. (C5) 

k 

We see that wi is the left eigenvector of WjNj with 
eigenvalue is the largest eigenvalue of wjNj, 

since the eigenvector has positive elements. As a result 
is common left eigenvector of Nj for all j’s, with eigen¬ 
value Wj . Since is non-negative, Wj is the largest eigen¬ 
value of Nj. Therefore 

Wj = dj (C6) 

is the quantum dimension of type-j particle (see 

eqn. (A40)). The largest left eigenvalue of A^i is 1. Thus 

di = Wi = 1 . 

We can represent a type-1 string by a dash line. By 
examine the 0-move with k = 1: 



we see that we can remove (or replace) any vertex with 
dash line by including a factor We can add a 

vertex with dash line by including a factor In other 

words, a vertex with dash line is just a weighting factor 
^- 1 / 2 ^ The unitary m-fusion categories with the trivial 
particle type will be called unitary fusion categories. 
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di 


FIG. 9: A loop of type-i string has an amplitude di. 


2. Amplitudes for loops 

The tensors characterize the four types 

of linear relations between graphs with some local differ¬ 
ences. Those local changes are almost complete, in the 
sense that any graphs of strings can be reduce to graphs 
that contain only isolated loops. Since the amplitude 
of a graph that contain disconnect parts is given by the 
product of the amplitudes for those parts, therefore, if we 
know the amplitude for single loops of string, then the 
amplitude of any string configuration can be computed 
from the tensor data 

With the presence of trivial particle type in the unitary 
fusion category, we can determine the amplitude for a 
loop of f-string. Using the rule of adding dash lines (the 
trivial strings) and O-move eqn. (A33), we find 



FIG. 10: A Y-move can change a linking of three loops in 
(a) to two linkings of two loops in (b). The amplitude for two 
linked i-loop and /-loop in (c) is denoted as S'l;"*'. 


The commutative unitary fusion category for planar 
graphs plus the extra structure for non-planar graphs and 
their amplitudes will give us a modular tensor category 
theory. In this section, we will derive many conditions 
that involve amplitudes of non-planar graphs. 



We may choose A = 1. This allows we to deter¬ 

mine (see Fig. 9) 

A (O) = 

Appendix D: Modular tensor category for the 
amplitudes of non-planar string configurations 

1. Commutative unitary fusion category 


2. Amplitude for linked loops and Verlinde formula 

As the first application of non-planar graphs, consider 
a three linked loops in Fig. 10a. We can evaluate the 
graph in two ways: (a) we fuse the i-loop and j-loop 
using eqn. (A4) to produce a single fc-loop; (b) we use 
a Y-move to change the linking of three loops to two 
linkings of two loops. We defined the amplitude for two 
linked loops i and j as (see Fig. 7). satisfies 

(D3) 

This allows us to obtain 

cLnk oLnk 

E , (D4) 

k * 

which is the tensor category version of Verlinde formula. 


We have being considering planar graphs and the re¬ 
lated fusion category theory. In this section we will con¬ 
sider non-planar graphs. Since the particles now live in 
2D space, the fusion of the particles satisfies 


j = j 

(DI) 

K' = 

(D2) 


So the fusion of 2D particles are commutative (while the 
fusion of ID particles may not be commutative). The 
fusion with Nl^ = is called commutative. Also, we 
assume the existence of trivial particle type. Thus, in 
this section, the fusion of the particles is described by a 
commutative unitary fusion category. 


3. Degenerate ground states on torus and 
excitation basis 

To obtain the algebraic structure for non-planar 
graphs, let us first try to represent the degenerate ground 
states on torus graphically. One of the degenerate ground 
state that corresponds to the trivial quasiparticle i = 1 
can be represented by an empty solid torus x Dyf. (see 
Fig. 11a), where the circle in y-direction Sy is a bound¬ 
ary of the disk Dy^. In other words, the path integral on 
the space-time 5”,), x Dy^. give rise to the state |1) on the 
surface Si x S^. We denote such a state as |1). Other 
degenerated ground states can be obtained by the action 
of the Wf operators 

|*) = VFf|l). (D5) 
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FIG. 11: (a): The ground state |1) on a torus that corre¬ 

sponds to the trivial quasiparticle can be represented by an 
empty solid torus, (b): The other ground state |i) that cor¬ 
responds to a type i quasiparticle can be represented by an 
solid torus with a loop of type i in the center. 



FIG. 12: The graphic representation of Wf\i). 


|i)’s form a orthonormal basis if we assume {WfY = Wf 
and (i|l) = Sn. This is because 

{j\i) = (i| = {i\wjw^\i) = 

k 

= (l|iVf |1) = d,, (D6) 

We will call such a basis of the degenerate ground state 
an excitation basis. 

Since |1) is created by the tunneling operator , |z) 
can be represented by adding a i-loop that corresponds 
to the Wf operator to the center of the solid torus (see 
Fig. lib). 

|i) is a natural basis, where the matrix elements of 
and Wf have simple forms. From eqn. (A59), we see that 

WfWf\l) = Wflz) = |fc). (D7) 

k 

The action of VFJ on |*) is represented by Fig. 12. From 
Fig. 13, we find that 

yrii ^Lnk 

WJ\i) = = ^\i). (D8) 






FIG. 14: Gluing two solid tori D^t x Sy without twist forms 
a. X . The gluing is done by identifying the {x,y) point 
on the surface of the first torus with the (a;, —y) point on the 
surface of the second torus. If we add an additional S twist, 
i.e. if we identify {x, y) with {—y, —x), the gluing will produce 
a S\ 



FIG. 15: (Golor online) A “self-loop” with canonical framing 
corresponds to a twist by 2?!. A twist by 27r induces a phase 
ei27rsi dehnes the spin Si of the particle. 


We see that \i)’s are common eigenstates of the com¬ 
muting set of operators Wf. The corresponding eigen- 

gLnk 

value for Wf is -ff—. We see that different |l)’s have 
different set of eigenvalues, which support our assump¬ 
tion (i|l) = Sii. 


4. The relation between and S 

The amplitude of two linked loops, and the rep¬ 

resentation of the modular group in the excitations basis, 
S, are closely related. From Fig. 14, we see that 

Sij = (*|5|j) = (two linked loops in S^) = 

(D9) 

where is the volume independent part of the par¬ 
tition function on three-sphere Since S is unitary 

and = di, we see that 

ZinAS^) = l/D. (DIO) 



FIG. 13: (a): The graphic representation of Wj\i). (b): 

The graphic representation of Wf\i). (c): The Y-move can 
deformed the graph in (b) to the graph in (c). The shaded 
area represents the hole of the torus. 


FIG. 16: The two “self-loops” in (a) are “right-handed” and 
correspond to the same twist. The two “self-loops” in (b) are 
“left-handed” and also correspond to the same twist that is 
opposite to that in (a). 
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TABLE XI: A list of primitive bosonic topological orders in 2+lD with rank N — 2,3, - ■ ■ ,6. The list contain all topological 
orders with rank N = 2,3, 4. The list may not be complete for rank N = 5,6. However, it contains all rank N — 5 primitive 
topological orders with < 120, and all rank N — 6 primitive topological orders with < 60. 


iVf 

5'top 


di,d2, ■ ■ ■ 

■Sl, S2, • • • 

Jf-matrix/SCA 

2f 

0.5 

2 

1,1 



0,i 

(2) 

2^1 

0.5 

2 

1,1 



0,-i 

(-2) 

^14/5 

0.9276 

3.6180 

l,Cl 



o.i 

(Ai, 3 ) 1 / 2 , (G 2 ,1) 

9^ 

^-14/5 

0.9276 

3.6180 

1,C3 



0 —- 
U, 5 

(Ai,— 3 ) 1 / 2 , (G 2 ,—1) 

3f 

0.7924 

3 

1,1,1 



oil 

3 ’ 3 

(2 2;1) 

3^2 

0.7924 

3 

1,1,1 



0 -1 -1 

-(2 2;1) 

oB 

^Z/2 

1 

4 

1,1,^ 



0 1 A 

2 ’ 16 

(Ai,2),(B9,1) 

oB 

'^5/2 

1 

4 

1,1,^ 



0 1 A 

2 ’ 16 


oB 

'^7/2 

1 

4 

1,1,^ 




(B3,1) 

oB 

'J-7/2 

1 

4 

1,1,^ 



2 ’ 16 

(B4,1) 

oS 

^-5/2 

1 

4 

1,1,^ 



0 1 -A 

(B5,1) 

oB 

'J-3/2 

1 

4 

1,1,^ 



0 1 -A 

(B6,1) 

oB 

^-1/2 

1 

4 

1,1,^ 



0 1 -A 

57,1) 

oB 

^1/2 

1 

4 

1,1,^ 



0 ^ ^ 

(58,1) 

CO 

“4 

1.6082 

9.2958 

1,C5,C5 



0-12 

(Ai, 5 ) 1/2 

oB 

’^-8/7 

1.6082 

9.2958 

i,Ci,Cl 



0 1 -2 
’ 7 ’ 7 

(Ai, — 5 ) 1/2 


1 

4 

1 , 1 , 1,1 



0 , 0 , 0,1 

(0,0; 2) 

4f 

1 

4 

1 , 1 , 1,1 



0 1 i i 

(4) 

4?i 

1 

4 

1 , 1 , 1,1 



0-1-11 

8 ’ 8 ’ 2 

(-4) 

4f 

1 

4 

1 , 1 , 1,1 



0 13 1 

(2 2 2;1 1;1) 

4^3 

1 

4 

1 , 1 , 1,1 



0 -3 _3 1 

8 ’ 8 ’ 2 

-(2 2 2;1 1;1) 

4f 

1 

4 

1 , 1 , 1,1 



0 , 1 5 

(2 2 2 2;1 0 0;1 0;1) 

^10/3 

2.1328 

19.234 

1,C7,C?,C7 



0 1 2 _1 

(Ai, 7 ) 1/2 

A^ 

^-10/3 

2.1328 

19.234 

l,C7\C7",Cl 



O’-l’-f’ 5 

(Ai,-7)i/2,(G2,2) 

5® 

1.1609 

5 

1,1, 1,1,1 



oil 1 1 

5 ’ 5 ’ 5 ’ 5 

(2 2; 3) 

5f 

1.1609 

5 

1,1, 1,1,1 



r, 2 2 2 2 

5 ’ 5 ’ 5 ’ 5 

(2 2 2 2;1 1 0;1 0;1) 

5?’“ 

1.7924 

12 

1,1, V3, ^,2 



0 0 i i 

8 ’ 8 ’ 3 

(Ai,4) 

5?’" 

1.7924 

12 

1,1,V3, ^,2 



8 ’ 8 ’ 3 

(5f’“KI2f KI2 ®i)i/4 

5^2“ 

1.7924 

12 

1,1,V3, ^,2 



00 1 ^ 1 

(Ai,-4) 


1.7924 

12 

1,1,V3, ^,2 




(5(?^“ Kl 2f Kl 2 ®i)i/4 

^16/11 

2.5573 

34.646 

1 , C 9 , cl, cl, cl 



r, 2 2 1 5 

11 ’ 11 ’ 11 ’ 11 

(Ai, 9 ) 1 / 2 , ( 54 , 2 ) 

^-16/11 

2.5573 

34.646 

1 , cl, cl, cl, cl 



n 2 2 15 

11 ’ 11 ’ 11 ’ 11 

(Ai,— 9 ) 1 / 2 , (58,3) 

"^18/7 

2.5716 

35.342 

l,C5,C5,Cl2,Cl2 



0 1 113 

(Ai, 12 ) 1 / 4 , (A 2 , 4 ) 1/3 

r:B 

'^-18/7 

2.5716 

35.342 

1, Cs, cl, C 12 , C 12 



0 1 1 _i -3 

7 ’ 7 ’ 7 ’ 7 

(A 3 , 3 ) 1/4 

nB ,a 

6 o 

2.1609 

20 

1 ,1, 2, 2, V5,V5 



0 0 1 —1 0 2 

( 55 , 2 ) 1 / 4 , ([ 7 ( 1 ) 5 / 12 ) 1/2 

CiB,b 

6 o 

2.1609 

20 

1 ,1, 2, 2, ^/5,^/5 



0 , 0 , 5 , 5 > 4 > 4 

(6^’“KI2f KI2 (!i)i/4 

nB ,a 

64 

2.1609 

20 

1,1,2,2,V5,V5 



002 -21-1 

(52,2) 

aB,b 

64 

2.1609 

20 

1 ,1, 2, 2, V5,V5 



0 0 2 —2 0 1 

VJ, u, 5 , 5 ? 2 

(6fKl 2f Kl 2 (!i)i/4 

^46/13 

2.9132 

56.746 

1, Cii, Cii, Cii, Cii, Cii 



Q 4 2 6 6 1 

13 ’ 13 ’ 13 ’ 13 ’ 13 

(Ai, 11 ) 1/2 

^-46/13 

2.9132 

56.746 

1, Cii, Cii, Cii, Cii, Cii 



0 4 2 6 6 1 

13 ’ 13 ’ 13 ’ 13 ’ 13 

(Ai, — 11 ) 1/2 

'^8/3 

3.1107 

74.617 

l,C7",C7",C7",Cl6,Cm 



0 1 1 1 1 _1 
’ 9 ’ 9 ’ 9 ’ 3 ’ 3 

(Ai, 16 )i/4 

6 - 8/3 

6 f 

6^2 

3.1107 

3.3263 

3.3263 

74.617 

100.61 

100.61 

l,C7^C7^C7^Cl6,Cf6 

1 3 +T 2 T 3+Tm 3+Vm 

-*-5 2 ’ 2 ’ 2 ’ 

1 3+V^ 3+V^ 3+V^ 

2 ’ 2 ’ 2 ’ 

5 +V 2 T 

2 ’ 
5+V^ 

2 ’ 

7+V^ 

2 

7+V^ 

2 

0 1 1 1 1 1 

9’ 9’ 9’ 3’ 3 
0 ,- 2 ,- 2 , |, 0 , 1 

0 1 2 _3 Q _1 

^>7175 7 ? 3 

(Ai,—16 )i/4, (A 2 , 6 ) 1/9 
(G2,-3) 

(G2,3) 
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TABLE XII: A list of primitive bosonic topological orders in 2+lD with rank N = 7,8,9. The list may not be complete. 
However, it contains all rank N = 7 primitive topological orders with < 30, all rank N = 8 primitive topological orders 
with < 20, and all rank N = 9 primitive topological orders with < 15. 



5'top 


di,d2, ■ ■ ■ 

Si, S2,■ • ■ 

TC-matrix/SCA 

rjB ,a 
' 2 

1.4036 

7 

1,1,1,1,1,1,1 

o 1 1 2 2 3 3 

7’ 7 

(4 4; 3) 

’’jB ,0. 

^ -2 

1.4036 

7 

1,1,1,1,1,1,1 

n 1 1 2 2 3 3 

7? 7? 7? 7’7’7 

-(4 4; 3), (A6,1) 

7B 

'9/4 

2.3857 

27.313 

Ij 1) Cb) Ce I cl 1 cl 1 cl 

PI 1 3 3 1 1 15 

2 ’ 32 ’ 32 5 4 ’ 4 ’ 32 

(4li,6) 

jB 

' 13/4 

2.3857 

27.313 

ij 1 ) cl) cl I cli cli cl 

0 17 7 1 1 13 

2 ’ 32 ’ 32 ’ 4 ’ 4 ’ 32 

(7g/4H4f)i/4 

jB 

'-15/4 

2.3857 

27.313 

1 ) 1 ) cl) cl) cl: cl: cl 

Q 1 11 11 1 1 9 

2 ’ 32 ’ 32 ’ 4 ’ 4 ’ 32 

(7/3/4 H 4f )i/4 

7B 

'-11/4 

2.3857 

27.313 

1 ) 1 ) cl) cl) cl: cl: cl 

1^ 1 15 15 1 1 5 

2 ’ 32 ’ 32 ’ 4 ’ 4 ’ 32 

(7/45/4 H4f)i/4 

7B 

^ -7/4 

2.3857 

27.313 

1 ) 1 ) cl) cl) cl) cl) cl 

PI 1 13 13 1 1 1 

(7/41/4 H4f)i/4 

7B 

^ -3/4 

2.3857 

27.313 

1 ) 1 ) cl) cl) cl) cl) cl 

0 1 _ 9 _ 9 1 _1 3 

2 ’ 32 ’ 32 ’ 4 ’ 4 ’ 32 

(7/,/4H4f)i/4 

7B 
' 1/4 

2.3857 

27.313 

1 ) 1 ) cl) cl) cl: cl: cl 

o 1 5 5 1 17 

2 ’ 32 ’ 32 ’ 4 ’ 4 ’ 32 

(7/3/4 H4f)i/4 

7B 
' 5/4 

2.3857 

27.313 

1 ) 1 ) cl) cl) cl: cl: cl 

1 1 11 1 11 

2 ’ 32 ’ 32 ’ 4 ’ 4 ’ 32 

(7f/4H4f)i/4 

7B 

' 7/4 

2.3857 

27.313 

1 ) 1 ) cl) cl) cl) cl) cl 

PI 1 13 13 1 11 

2 ’ 32 ’ 32 5 4 ’ 4 ’ 32 

(Ce, 1) 

7B 
• 11/4 

2.3857 

27.313 

1 ) 1 ) cl) cl) cl) cl) cl 

PI 1 15 15 1 15 

2 ’ 32 ’ 32 ’ 4 ’ 4 ’ 32 

(7//4H4f)i/4 

7B 
• 15/4 

2.3857 

27.313 

1 ) 1 ) cl) cl) cl: cl: cl 

1 11 11 1 19 

2 ’ 32 ’ 32 ’ 4 ’ 4 ’ 32 

(7/1/48 4f)i/4 

7B 

'-13/4 

2.3857 

27.313 

1 ) 1 ) cl) cl) cl: cl: cl 

n 1 7 7 1 1 13 

2 ’ 32 ’ 32 ’ 4 ’ 4 ’ 32 

(7/5/4 8 4/)i/4 

7B 

^ -9/4 

2.3857 

27.313 

1 ) 1 ) cl) cl) cl) cl) cl 

PI 1 3 3 1 1 15 

('^-13/4 ® di )l/4 

7B 

^-5/4 

2.3857 

27.313 

1 ) 1 ) cl) cl) cl) cl) cl 

PI 1 1 1 1 1 11 

2 ’ 32 ’ 32 ’ 4 ’ 4 ’ 32 

(7/g/4 8 4/)l/4 

7B 

'-1/4 

2.3857 

27.313 

1 ) 1 ) cl) cl) cl) cl) cl 

0 15 5 1 1 7 

2 ’ 32 ’ 32 ’ 4 ’ 4 ’ 32 

(7/5/4 8 4f)i/4 

7B 

'3/4 

2.3857 

27.313 

1 ) 1 ) cl) cl) cl: cl: cl 

0 19 9 1 1 3 

2 ’ 32 ’ 32 ’ 4 ’ 4 ’ 32 

(7/4/4 8 4/)i/4 

rjB,b 

'2 

2.4036 

28 

1,1,2,2,2,77,77 

7 ? 7 1 7 ’ 8 ’ 8 

{U{l)r/Z2)i/2 

fjB ,c 

' 2 

2.4036 

28 

1,1,2,2,2,77,77 

0 0 1 2 3 1 3 

^^5 7 ? 7 ? 7 ? 8 ? 8 

(7/’^ 8 2/82/1)1/4 

rjB,b 
' -2 

2.4036 

28 

1,1,2,2,2,77,77 

0 0 1 2 3 13 

7 ? 7 ’ 7 ’ 8 ’ 8 

(B3,2), (D7, 2)1/2 

i-wB ,c 
' -2 

2.4036 

28 

1,1,2,2,2,77,77 

0 0 ^ 2 3 1 3 

(7//8 2/82/1)1/4 

7B 

'8/5 

3.2194 

86.750 

1) Cl3) Cl3: Cl3) Cl3: Cl3) Cl3 

0 _1 A 0 2 1 _1 

5 ’ 15 ’ 5’ 3’ 5 

(Ai, 13)1/2 

7-8/5 

3.2194 

86.750 

1) Cl3: Cl3) Cl3: Cl3) Cl3) Cl3 

0 1 -A 0 -2 -i 1 

(Ai,-13)1/2 

7f 

3.2715 

93.254 

1) cl : cl ) 1 + cl : f + C| ) 2C|: 1 + 2C| 

Q 1 1 1 1 _3 0 

^’2’2’4’4’ 8’ 

(A2,5)1/3 

7?i 

3.2715 

93.254 

1,C|)C|,1 + C|,1 + C|,2C|,1 + 2C| 

0) |, “i: |) 0 

(A2, —5)1/3 

8f’“ 

1.5 

8 

1)1) 1)1,1)1) 1)1 

0 0 1 1 1 1 7 7 

16 ’ 16 ’ 4 ’ 4 ’ 16 ’ 16 

(8) 

8/’" 

1.5 

8 

1)1) 1)1,1)1) 1)1 

o Q 3 3 115 5 

16 ’ 16 ’ 4 ’ 4 ’ 16 ’ 16 

(2 2 2;1 2; 3) 

8fl“ 

1.5 

8 

1)1) 1)1,1)1) 1)1 

16 ’ 16 ’ 4 ’ 4 ’ 16 ’ 16 

(-8) 

8®/ 

1.5 

8 

1,1,1,1,1,1,1,1 

Q Q 3 3 1 1 5 5 

16 ’ 16 ’ 4 ’ 4 ’ 16 ’ 16 

-(2 2 2;1 2; 3) 

QB,a 

2.5849 

36 

1,1, 2, 2, 2, 2, 3, 3 

000-—- - 0- 

(S4,2) 

oB,b 

2.5849 

36 

1,1, 2, 2, 2, 2, 3, 3 

0 0 0 ^2 4 Q 1 

(B4,-2) 

^B,c 

2.5849 

36 

1,1, 2, 2, 2, 2, 3, 3 

0 0 0 ^ 241 1 

9 ’ 9 ’ 9 ’ 4 ’ 4 

(8/’“ 8 2/82/1)1/4 

qB ,d 

2.5849 

36 

1,1, 2, 2, 2, 2, 3, 3 

0 0 0 12 4 1 1 

(8®’^ 8 2/82/1)1/4 

®62/17 

3.4879 

125.87 

1) ClS) Cl5) Cl5) Cl5: Cl5) Cl5) Cl5 

0 5 2 8 6 4 5 3 

17’ 17’ 17’ 17’ 17’ 17’ 17 

(Ai, 15)1/2 

®-62/17 

3.4879 

125.87 

1) Cl5) Cl5) Cl5) Cl5: Cl5) Cl5) Cl5 

n5 2 8 645 3 

17’ 17’ 17’ 17’ 17’ 17’ 17 

(Ai,-15)1/2 

^0 

1.5849 

9 

1,1, 1,1, 1,1, 1,1,1 

000-------- 

(4 4; 5) 

^0 

1.5849 

9 

1,1, 1,1, 1,1, 1,1,1 

0 0 0 1 1 2 2 4 4 

9 ’ 9 ’ 9 ’ 9 ’ 9 ’ 9 

-(4 4; 5) 



This way, we ontain an important relation: 


FIG. 17: A figure “8” of type-i string has an amplitude 




(Dll) 


that connects the amplitude of two linked loops to a mod¬ 
ular transformation of the degenerate ground state on 
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FIG. 18: Gluing two solid tori x with an additional 
T twist, i.e. identifying {x,y) with (x + y,—y), will produce 
a 5"^ X . The i-loop in y-direction in the second solid torus 
at right can be deformed into a i-loop in the first solid torus 
at left. We see that the loop is twisted by 2n in the anti¬ 
clockwise direction. 


(a) (b) 



(c) 


FIG. 19: A double figure “8” in (a) is equal to two linked 
figure “8” in (b) after sliding one of the figure “8”. Applying 
two Y-moves to (b), we obtain (c) 

torus. This allows us to rewrite eqn. (D4) as 

= (D12) 

k 

which is the Verlinde formula. 


5. The spin Si of topological excitations and T 

We have studied two linked loops, which is a string 
configuration with crossing. Another important string 
configuration with crossing is a “self-loop” (see Fig. 15). 
Such a “self-loop” corresponds to a twist by 27r, which 
equal to a straight line with a phase . Here Si is the 
spin of the type-i topological excitation, which is defined 
mod 1. We also note that the handness of the “self-loop” 
determines the direction of the twist (see Fig. 16). As 
a result, a figure “8” of type-f string has an amplitude 

e27riSi^. pjg_ 27). 

It is clear that the Dehn twist T, when acting on 
|z) = Wf\l), will twist the string i by 27r and induces a 
phase (see Fig. 18). The Dehn twist T also change 



the space-time metrics which may causes an additional 
i-independent phase, which is denoted as This 

way, we ontain another important relation: 

T|z) = e-2"*^e2"‘"*|f). (D13) 

6. Relation between and Sij 

To understand the relation (33), let us compute the 
amplitude of a double figure “8” in two ways, as shown 
ih Figs. 19 and 20. This allows us to show 

= (D14) 

a,k 

which can be simplified to 

(D15) 

dj di ^ 

or 

(D16) 

k 

Using S*-. = Sij, the above becomes eqn. (33). 

In this paper, we have derived most of the ,Si,c) 

conditions, expect the condition eqn. (28). Here, we 
would like mention that the condition eqn. (28) can be 
found in Ref. 104. 


Appendix E: List of primitive topological orders 

In this section, we give lists that contain more topolog¬ 
ical orders (see Tables XI and XH, where only the primi¬ 
tive topological orders are listed). The lists are generated 
using a different numerical code. 

The abelian states with di = 1 are de¬ 

scribed by AT-matrices. We use the notation 
{KiiK 22 -■ ■ ■,Ki2K23-■ ■ ■,Ki3K24:-■ ■ ]■ ■ ■) to denote 
the AT-matrices. 

In Ref. 106, we show that the non-abelian states can 
be generated by simple current algebra (SCA)^"^’®^’®^. 
(See also https://www.math.ksu.edu/~gerald/voas/ ) 
The SCA’s are denoted by {R,±k)a (see Ref. 103,106), 
where R = An, Bn,Cn, D^, etc , and {R,—k)a is the 
time-reversal conjugate of {R, +k)a. The last column of 
Tables XI and XH indicates how the corresponding topo¬ 
logical order is realized by the AT-matrix state or the SCA 
state. 




FIG. 20: A double figure “8” in (a) is changed to (b) after a 
Y-move. (b) is equal to (c) after sliding the vertices. 
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